METHOD OF ANALYZING ELECTROMAGNETIC INTERFERENCE 



Background of the Invention 

1. Field of the Invention 

The present invention relates to a method of 
analyzing electromagnetic interference (EMI) , and more 
particularly, to a method of analyzing electromagnetic 
interference arising in a large-scale, high-speed LSI 
(large-scale integrated circuit) by means of high-speed, 
highly-accurate logic simulation. 

2 . Desaription of the related Art 

As is often the case, in the field of computers, LSIs 
find a broadening range of applications, from 
communications devices, such as cellular phones, to 
general household products, toys, and automobiles. 
Electromagnetic interference arising in such a product 
induces radio interference noise to arise in a receiver, 
such as a TV set or a radio, or faulty operations of another 
system. In order to prevent this problem, the entirety 
of a product is shielded, or filters are provided in a 
product. With a view towards preventing an increase in 
the number of components and cost and difficulty in 
preventing occurrence of electromagnetic interference in 
a product, strong demand exists for suppression of noise 
in an LSI alone. 
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Under such a situation, an LSI is ranked as a key 
device for any product which contains an LSI . Demand 
exists for an larger-scale, high-speed LSI for ensuring 
competitiveness of a product. In a situation in which the 
5 cycle of product development becomes shorter, design- 
automation of an LSI is indispensable for satisfying the 
demand. There is growing necessity for adopting 
synchronous circuit as a condition for introducing a 
state-of-the-art design-automation method. In a case 
10 where all circuits of a large-scale, high-speed LSI 
operate synchronously with a reference clock signal , 
instantaneously-changing current becomes very large, 
which induces an increase in electromagnetic 
interference . 

15 The present invention relates to a simulation method 

which enables evaluation of EMI indispensable for reducing 
electromagnetic interference while maintaining a tendency 
toward a larger-scale, higher-speed LSI . 

Noise imposed on another device by an LSI is roughly 

20 classified into two types; radiation noise, and conduction 
noise. Radiation noise emanated directly from an LSI 
includes noise emitted from internal wires of an LSI. 
However, internal wires do not act as an antenna of large 
size. As a matter of course, it is considered that the 

25 noise emitted directly from an LSI will pose a problem in 
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the future, in association with an improvement in the 
operation frequency of an LSI. However, as of now, the 
noise emitted from the inside of an LSI is considered 
trivial - 

5 In contrast, conduction noise affects another device 

mounted on a printed board by way of direct 
interconnections, such as internal wires of an LSI, lead 
frames, or wiring patterns provided on a printed wiring 
board. Noise is emitted from such interconnections while 

10 the interconnections are taken as the source of 

origination or as an antenna. The antenna constituted of 
the interconnections is much larger than that constituted 
by internal wires of an LSI and is a dominant element in 
terms of electromagnetic emission. 

15 A power line and a signal line can act as paths along 

which conduction noise developing in an LSI travels. In 
consideration of an electromagnetic field in the vicinity 
of an LSI, noise which results from variation in an 
electric current of a power source being emitted from a 

20 power line serving as an antenna is considered to be 

dominant. There may be a case where a ringing overshoot 
stemming from variation in a signal poses a problem. 
However, there more frequently arises a case where 
variation in an internal power level of an LSI propagates 

25 as a signal waveform, to thereby presents a problem. Noise 
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emitted from a power line or a signal line is considered 
to have a strong correlation with variation in the electric 
current of a power source (hereinafter referred to as a 
* source current"). 
5 A source current of a CMOS circuit will now be 

described by reference to a simple inverter circuit. In 
a case where variation arises in a voltage applied to an 
inverter circuit, there flows a load capacitance 
charge/discharge current, which is the primary source 

10 current of the CMOS circuit. In addition, a short circuit 
current flows together with the load capacitance 
charge/discharge current. In design of such a CMOS 
circuit, all circuits of an LSI are synchronized in 
accordance with constraints on the use of a design- 

15 automation tool. As a result of all circuits being 
synchronized, all circuits of the LSI operate 
simultaneously, and a peak current arises in a power source 
in synchronism with a reference clock signal. Further, 
in order to increase operating speed, or shorten a cycle, 

20 of the LSI, the capacitance of a transistor is increased 
so as to enable a charging/discharging operation to be 
completed within a short period of time. Eventually, a 
peak current increases. As a matter of course, the total 
source current of an LSI is increased even when the scale 

25 of an LSI is increased. Thus , the peak current of the power 
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source is increased, thereby inducing occurrence of an 
abrupt change in a source current. Such an abrupt change 
induces an increase in higher harmonic components, thereby 
resulting in an increase in electromagnetic interference. 
5 Highly-precise simulation of change in a source 

current, which may be said to primarily account for 
electromagnetic interference, is considered to be 
effective in evaluation of electromagnetic interference 
arising in an LSI. 
J7; 10 A current simulation method for effecting 

zl transistor-level current analysis, as will be described 

I** below, has conventionally been employed. 

Cl FIG. 116 is a block diagram showing the flow of 

CI processing operations pertaining to a conventional EMI 

g 15 analysis method. According to this method, on the basis 

of layout information pertaining to an LSI to be analyzed 
through use of a method of analyzing an electric current 
of a transistor, there is performed processing pertaining 
to a layout parameter extraction (hereinafter referred to 
20 simply as an "LPE") step 03 . Subsequently, there is 

performed processing pertaining to a circuit simulation 
step 06 for effecting simulation of a circuit on the basis 
of a switch-level netlist; a source-of -current modeling 
step 08; a source line LPE step O10; a transient analysis 
25 simulation step 012 ; and a FFT step 014 . 



5 




Processing pertaining to each of the foregoing 
processing steps will now be described by reference to FIG . 
116. 

First, there are input data Ol pertaining to the 
5 layout of a semiconductor integrated circuit (hereinafter 
referred to simply as "layout data Ol") , which circuit is 
to be subjected to EMI analysis; parameters of elements, 
such as transistor elements or various parasitic wiring 
elements (e.g. , resistors and capacitors) ; and an LPE rule 

10 02 for defining a form in which extracted layout parameters 
are to be output. In step 03 , parameters of the respective 
elements included in the layout data Ol are calculated in 
accordance with the LPE rule 02 , to thereby prepare a 
netlist 04 . In step 03 , parameters of parasitic elements 

15 of a power source (and the ground) are not extracted. 

In step 06 are input the netlist 04 prepared in step 
03 and a test pattern OS for replicating a desired logic 
operation in a circuit serving as an object of analysis. 
There are calculated a load capacitance charge/discharge 

20 current and a short circuit current, both responding to 
the operating state of an internal circuit, wherewith 
information 07 about the waveform of an electric current 
of each transistor (hereinafter referred to simply as 
"current waveform information 07") is produced. The 

25 first operation of the processing pertaining to step 06 
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is effected on the assumption that the potential of a power 
source (and that of ground) is a variation-free, ideal 
potential . 

In step 08 , the current waveform information pieces 
5 07 concerning the respective transistors prepared in step 
06 are entered. Each of the thus-entered current waveform 
information 07 pieces is modeled into a mode which can be 
applied to subsequent step 012, wherewith information 09 
concerning a source-of -current element model (hereinafter 

10 referred to simply as " source-of -current element model 
information 09") is prepared. According to a common 
method, in order to alleviate the load associated with 
processing pertaining to step 012 and subsequent steps, 
each functional circuit block consisting of a plurality 

15 of transistors is modeled as a source-of -current element. 

Processing pertaining to step O10 differs from 
processing pertaining to step 03 , only in that parameters 
of parasitic elements of a power source and parameters of 
those of parasitic elements of ground wiring (e.g., 

20 resistors, decoupling capacitance , and 1 ike elements ) are 
extracted and subjected to EMI analysis in lieu of 
extracted parameters of transistor elements and extracted 
parameters of various parasitic wiring elements. Hence, 
repeated explanation of processing pertaining to step O10 

25 is omitted. In step O10, a power source (and ground) 
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wiring netlist Oil is prepared. 

In step 012 , there are entered the source-of -current 
element model information 09 prepared in step 08 , the power 
source (and ground) wiring netlist Oil prepared in step 
5 OlO, and an impedance 016 (including, resistance, 

capacitance, and inductance) of a wire or lead frame. 
Through analysis of these input data carried out by a 
transient analysis simulator typified by SPICE , variation 
in the line voltage of a circuit to be analyzed is 



prepared . 

Subsequently, processing pertaining to step 06 is 
performed again. In contrast with the first operation of 
the processing pertaining to step 06 having been effected 

15 on the assumption that the potential of the power source 
(and the ground) is a variation-free, ideal potential, the 
line voltage drop result 016 prepared in step 012 is 
entered, and the current waveform information pieces 07 
concerning the respective transistors are prepared again 

20 in consideration of the line voltage drop. Similarly, 
processing pertaining to steps 08 and 012 is performed 
again . 



effected several times in a looped manner, wherewith there 
25 is produced a current waveform result 013 which duplicates 




Processing pertaining to steps 06, 08 , and 012 is 
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variation in a line voltage with high precision. 

In step 014 , the current waveform result 013 prepared 
in step 012 is entered and subjected to fast Fourier 
transformation (hereinafter abbreviated FFT) , to thereby 
5 enable frequency spectrum analysis. There is obtained an 
EMI analysis result 015. 

In the conventional example, the precision of 
verification varies greatly according to combination of 
processing pertaining to the LPE step 03 , processing 

10 pertaining to the source line LPE step O10, and processing 
pertaining to the source-of -current modeling step 08 . 
However, a certain level of accuracy of analysis can be 
expected. A transient analysis simulator typified by 
SPICE is used for transistor-level analysis of an electric 

15 current. Hence, a limitation is imposed on the level of 
a circuit to be analyzed, and an enormous amount of 
processing time is required. The level of a semiconductor 
integrated circuit has increased recently, and 
establishment of an EMI analysis method which enables 

20 high-speed analysis of an electric current on a larger 
scale than a transistor level by dealing with the analysis 
in abstract more than a transistor level dealing is 
desired . 

A gate-level current analysis method (hereinafter 
25 referred to as a "gate-level power consumption analysis 
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method") has conventionally been proposed as a current 
analysis method which can be made faster . This gate-level 
current analysis method is used for analyzing power 
consumption . 

5 The gate-level power consumption analysis method for 

estimating power consumption on a gate level comprises the 
following steps : namely, a step of preparing beforehand 
a library from the total amount of electric current 
required per change when a logic signal appearing at an 

10 output port of a circuit element is changed from a logic 
value of 0 to 1 or from a logic value of 1 to 0; and a step 
of adding, to a total amount of current of the entire 
circuit at each time, a total amount of electric current 
required per change when a logic signal appearing at an 

15 output port of an instance in a circuit, the circuit being 
an object of estimation of power consumption, is changed 
from a logic value of 0 to 1 or a logic value of 1 to 0, 
to thereby estimate the total power consumed by the overall 
circuit (i.e., a value obtained by means of multiplying 

20 the total amount of current by an external voltage) . FIG. 
109 shows the configuration of an example system for 
effecting the conventional, popular gate-level power 
consumption analysis method. The gate-level power 
consumption analysis system shown in the drawing is 

25 operated through use of an amount-of -current analyzer 
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comprising circuit connection information storage means 
1001; signal pattern storage means 1002; means 1003 for 
storing the amount of electric current flowing through an 
element (hereinafter referred to simply as "element 
5 current storage means 1003") ; to tal-amount-of -current 
storage means 1004; and total -amount-of -current 
calculation means 1005. 

Individual elements constituting the gate-level 
current analysis system shown in FIG. 109 will now be 

10 described. Procedures for calculating the amount of 

current shown in FIG. 113 will be described, by reference 
to the gate-level circuit connection information shown in 
FIG. 110, a pattern of signal change arising at each 
terminal shown in FIG. Ill, and the amount of current 

15 flowing through respective elements shown in FIG. 112. 

The circuit connection information storage means 
1001 stores circuit connection information concerning 
connection of circuits which are objects of analysis. 
Circuit connection information, such as that shown in FIG. 

20 110, is stored in the circuit connection information 

storage means 1001 in advance. The circuit connection 
information is generally called a netlist. 

The circuit connection information comprises an 
external input port D, an external clock signal input port 

25 CLK, an external output port Q, flip-flops FF1 and FF2 , 
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a buffer element BUF1 , and wires Dl and Ql for 
interconnecting elements. Each of the flip-flops FF1 and 
FF2 comprises a data input pin D, a clock signal input pin 
CK, and a data output pin Q. A logic value appearing at 
5 the data input pin D immediately before arrival of a 

leading edge of a clock signal input to the clock input 
pin CK is transmitted to the data output pin Q. 

The buffer element BUF1 has a data input pin A and 
a data output pin Y. A logic value appearing at the data 

10 input pin A is output to the data output pin Y. 

The above elements have a lag, such as an intra- 
element lag, and permit propagation of a signal with a time 
lag of 1 ns. The external input port D is connected 
directly to the input pin D of the flip-flop FF1 . The 

15 external clock input port CLK is a terminal for receiving 
variation in a clock signal from the outside and is 
connected directly to the input pin CK of the flip-flop 
FF1 and to the input pin CK of the flip-flop FF2 . The 
external output port Q is connected directly to the output 

20 pin Q of the flip-flop FF2 and outputs a signal value of 
the flip-flop FF2 . A wire Dl interconnects the output pin 
Y of the buffer BUF1 and the input pin D of the flip-flop 
FF2 . A wire Ql interconnects the output pin Q of the 
flip-flop FF1 and the input pin A of the buffer BUF1 . 

25 The signal pattern storage means 1002 stores a 
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pattern of signal change applied to information concerning 
connections of a circuit to be analyzed. As shown in FIG. 
Ill, a pattern of change in a signal input to the input 
port CLK and a pattern of change in a signal input to the 
5 input port D shown in FIG. 110 are stored in the signal 
pattern storage means 1002 in advance. The horizontal 
axis in FIG. 110 represents the time at which change arises 
in a signal pattern. Upon each change, the signal pattern 

5 is changed between a logic value of 0 and a logic value 

M= 10 of 1. 

HI The element current storage means 1003 stores 

CO 

FU information concerning the amount of electric current 

H flowing through an element; that is, information 

|H concerning the total amount of electric current flowing 

if 15 when single change arises in an external terminal of each 

element. Information concerning the amount of electric 
current flowing through an element, such as that shown in 
FIG. 112, is stored beforehand in the element current 
storage means 1003. Reference numeral 1301 designates 
20 the amount of electric current flowing through a flip- 
flop FF when change arises in a signal input to the clock 
signal input terminal CK of the flip-flop FF . Reference 
numeral 1302 designates the amount of electric current 
flowing through the flip-flop FF when change arises in a 
25 signal output from the output terminal of the flip-flop 
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FF. Reference numeral 1303 designates the amount of* 
electric current flowing through the buffer BUF when 
change arises in a signal output from the output terminal 
Y of the buffer BUF. 

The total -amount-of -current calculation means 1005 
is given a flowchart such as that shown in FIG. 114. 

In step 1501, the circuit connection information 
which is stored in the circuit connection information 
storage means 1001 and is shown in FIG. 110 is read. In 
step 1502, a signal pattern stored in the signal pattern 
storage means 1002 and shown in FIG. Ill is read. In step 
1503, information concerning the amount of electric 
current (hereinafter referred to simply as "element 
current information'') stored in the element current 
storage means 1003 and shown in FIG. 112 is read. 

In step 1504, the signal pattern shown in FIG. Ill 
is delivered to the circuit connection information shown 
in FIG. 110, to thereby enable propagation of a signal. 
When change arises in a signal appearing at the external 
terminal of a circuit element written in the element 
current information, the amount of electric current 
flowing through the external terminal is added to the 
information concerning the total amount of electric 
current at each time. A signal propagation method used 
herein is identical with a signal propagation method 
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employed in a timing simulator. More specifically, the 
method is to calculate, in time sequence, change arising 
in a real circuit in accordance with a given signal pattern 
by means of tracking circuits which are affected by the 
signal pattern. The signal propagation method employed 
in a timing simulator is generally embodied in the form 
of a plurality of commercial tools. Hence, detailed 
description of the signal propagation method is omitted. 
FIG. 112 shows change in a logic value of a signal appearing 
at each pin, which has been caused as a result of 
propagation of the signal pattern in each step. 

FIG. 112 shows change in a logic value in the same 
form as that in which the signal pattern is shown in FIG. 
111. After completion of propagation of the signal, in 
step 1505 the information concerning the total mount of 
electric current calculated at respective times shown in 
FIG. 113 is stored in the total-amount-of -current storage 
means 1004. The total -amount-of -current storage means 
1004 stores the information concerning the total amount 
of electric current detected at respective times according 
to the foregoing procedures , in the form shown in FIG. 113. 

In a case where information concerning the total 
amount of electric current is calculated in greater detail , 
a load dependence coefficient of a current which is 
dependent on an output load capacitance (or on a charge 
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current) and the total amount of electric current which 
is not dependent on an output load capacitance (or on a 
short circuit current) are stored in an element current 
storage section. Thereafter, information concerning the 
capacitance of each wiring (i.e., capacitance 
information) is read, the sum of the load-independent and 
the product of the capacitance information and the load 
dependence coefficient is taken as the total amount of 
electric current. As described above, according to the 
conventional gate-level power consumption calculation 
method shown in FIG. 110, the amount of electric current 
flowing upon a single change arising in a signal appearing 
at an input or output pin is taken as a unit or as flowing 
momentarily. In other words, the conventional method 
determines only the total amount of electric current. In 
terms of power consumption, the accuracy of the thus- 
calculated total amount of current is sufficient. 
However, EMI analysis requires information concerning 
chronological change in electric current, and the total 
amount of electric current is not sufficient in terms of 
EMI analysis. 

The conventional example using the transistor-level 
current analysis method can be expected to yield a certain 
level of accuracy. However, a transient analysis 
simulator typified by SPICE is used for transistor-level 



16 




current analysis. A limitation is imposed on the level 
of a circuit to be analyzed, and an enormous amount of 
processing time is required- The scale of a semiconductor 
integrated circuit has increased recently, and 
5 establishment of an EMI analysis method which enables 
high-speed analysis of an electric current at a larger 
level than a transistor level is desired. 

In the conventional example using the gate-level 

O 

y3 current analysis method, high-speed analysis of a current 

H* 10 is feasible. However, only the total amount of electric 

HI current is determined. The gate-level current analysis 

GO 

fU method is sufficient in terms of power consumption but in 

= sufficient in terms of EMI analysis. 

H; Summary of the Invention The present invention has 

« 15 been conceived to solve the drawbacks of the conventional 

methods and is aimed at providing a method for enabling 
valuation of electromagnetic interference in an LSI 
through simulation within a realistic period of time. 

To this end, an EMI analysis method according to the 
20 present invention is characterized by means of adding a 
function of highly precisely analyzing a change in power 
current to a high-speed gate-level simulation. 

The present invention provides an electromagnetic 
interference analysis method for analyzing the amount of 
25 electromagnetic interference arising in an LSI by means 
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of performing a logic simulation, the method comprising: 

an instantaneous current calculation step of 
calculating the amount of instantaneous electric current 
from event information, the information being produced 
when a change arises in a signal and including the instance 
name of each cell in which the change has arisen, the name 
of the signal, a time at which the change has arisen, and 
transition information ; 

a modeling step of modeling the instantaneous 
electric current according to a predetermined rule; and 

an FFT processing step of subjecting to fast Fourier 
processing (hereinafter referred to as "FFT 
processing" ) the information concerning a change in 
electric current, the information being calculated 
through a modeling step. Here, the amount of 
instantaneous current means the amount of electric current 
induced by an event. The amount of electric current 
flowing through a cell when a change arises in an output 
of the cell as an event is taken as the amount of 
instantaneous electric current. 

Preferably, the modeling step includes an averaging 
step of averaging the instantaneous current over a certain 
discrete width, and the FFT processing step includes a step 
of subjecting to FFT processing information concerning a 
change in current, the information being produced by the 
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averaging step. 

Preferably, the modeling step includes a rectangular 
waveform modeling step of modeling the instantaneous 
current as a rectangular waveform whose height is 
calculated for information for each event such that the 
area of the rectangular waveform becomes equal to the 
amount of electric current of each event, and the FFT 
processing step includes a step of subjecting to FFT 
processing information concerning a change in current, the 
information being calculated in the rectangular waveform 
modeling step. 

Preferably, the modeling step includes a 
geometrically-similar rectangular waveform modeling step 
of modeling the instantaneous current as a 
geometrically-similar rectangular waveform whose height 
and width are calculated such that the area of the 
rectangular waveform becomes equal to the amount of 
electric current of each event, and the FFT processing step 
includes a step of subjecting to FFT processing 
information concerning a change in current, the 
information being calculated in the geometrically-similar 
rectangular waveform modeling step. 

Preferably, the modeling step includes a rectangular 
waveform modeling step of calculating the instantaneous 
electric current for information for each event, and a step 
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of modeling the instantaneous current as a rectangular 
waveform through use of the amount of electric current and 
a table representing the relationship between the width 
and height of a rectangular waveform, to thereby subject 
5 to FFT processing the information concerning a change in 
electric current calculated in the rectangular waveform 
modeling step. 

Preferably, the modeling step includes a step of 
calculating the instantaneous electric current for 

10 information for each event, and a rectangular waveform 
modeling step of modeling the instantaneous current as a 
rectangular waveform through use of a slew in input 
waveform and a table representing the relationship between 
the width and height of a rectangular waveform, to thereby 

15 subject to FFT processing the information concerning a 
change in electric current calculated in the rectangular 
waveform modeling step. 

Preferably, the modeling step includes a step of 
calculating the instantaneous electric current for 

20 information for each event, and a rectangular waveform 
modeling step of modeling the instantaneous current as a 
rectangular waveform through use of an output load 
capacitance and a table representing the relationship 
between the width and height of a rectangular waveform, 

25 to thereby subject to FFT processing the information 



20 




concerning a change in electric current calculated in the 
rectangular waveform modeling step. 

Preferably, the modeling step includes a step of 
calculating a drop in voltage from the amount of electric 
5 current flowing in each cell and the resistance of a power 
and correcting the amount of instantaneous electric 
current of each cell for each event, on the basis of the 
relationship between the drop in voltage and the amount 
of instantaneous electric current. 

10 Preferably, the modeling step includes . a step of 

calculating a drop in voltage from the amount of electric 
current flowing in each cell, the resistance of a power 
line, and the capacitance of an on-chip capacitor, and 
correcting the amount of instantaneous electric current 

15 of each cell for each event, on the basis of the 

relationship between the drop in voltage and the amount 
of instantaneous electric current. 

Preferably, the modeling step includes a step of 
transiently-analyzing a power RC of each cell and a cell 

20 power source, accurately calculating a drop in voltage, 
and a correction step of correcting the amount of 
instantaneous electric current of each cell for each event, 
on the basis of the relationship between the drop in 
voltage and the amount of instantaneous electric current. 

25 Preferably, the modeling step includes a triangular 
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waveform modeling step of modeling the instantaneous 
current as a triangular waveform which has a given width 
and whose height is calculated for information for each 
event such that the amount of instantaneous electric 
5 current becomes equal to the area of the triangular 

waveform, and the FFT processing step includes a step of 
subjecting to FFT processing information concerning a 
change in current, the information being calculated in the 
triangular waveform modeling step. 

10 Preferably, the modeling step includes a multi- 

order-function waveform modeling step of modeling the 
instantaneous current as a multi-order-function waveform, 
and the FFT processing step includes a step of subjecting 
to FFT processing information concerning a change in 

15 current, the information being calculated in the 
multi-order-function waveform modeling step. 

Preferably, the modeling step includes a exponential 
function waveform modeling step of modeling the 
instantaneous current as an exponential-function waveform, 

20 and the FFT processing step includes a step of subjecting 
to FFT processing information concerning a change in 
current, the information being calculated in the 
exponential-function waveform modeling step. 

Preferably, the modeling step includes a step of 

25 modeling the amount of instantaneous electric current 
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while separating the same into a short circuit electric 
current component and a charge current component. 

Preferably, the modeling step includes a calculation 
step of* calculating the height of a rectangular waveform 
5 from a library in which peak currents of cells are 
characterized according to the type of cell, and a 
rectangular waveform modeling step of modeling the amount 
of instantaneous electric current as a rectangular 
waveform, _the peak current calculated in the calculation 

10 step being taken as the height of the rectangular waveform 
and the area of the triangular waveform being equal to the 
amount of electric current of each event, and the FFT 
processing step includes a step of subjecting to FFT 
processing information concerning a change in current, the 

15 information being calculated in the rectangular waveform 
modeling step. 

Preferably, the calculation step includes a step of 
calculating a peak current for each cell from information 
concerning a slew in the cell, by reference to a library 

20 in which the relationship between a slew in input waveform 
and a peak current is characterized in the form of a table 
according to the type of cell . 

Preferably, the calculation step includes a step of 
calculating a peak current for each cell from information 

25 concerning a load capacitance of a cell, by reference to 
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a library in which the relationship between a load 
capacitance and a peak current is characterized in the form 
of a table according to the type of cell. 

Preferably, the calculation step includes a step of 
5 setting a plurality of peak currents for a composite cell 
and calculating the heights of a plurality of rectangular 
waveforms through use of a characterized library, and the 
rectangular waveform modeling step corresponds to a step 
of modeling the amount of electric current into a plurality 
10 of rectangular waveforms. 

Preferably, the calculation step includes a step of 
setting a plurality of peak currents for each of the rise 
and fall of a flip-flop (FF) cell and calculating the 
heights of a plurality of rectangular waveforms through 
15 use of a characterized library, and the rectangular 

waveform modeling step corresponds to a step of modeling 
the amount of electric current into a plurality of 
rectangular waveforms. 

Preferably, the calculation step includes a step of 
20 calculating the height of a rectangular waveform through 
use of a library in which peak currents are characterized, 
in consideration of the state of an input signal. 

Preferably, the modeling step includes 

a step of calculating a drop in voltage from the 
25 amount of electric current determined according to the 
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type of cell and from the resistance of a power line; and 
a correction step of characterizing, for each cell, 
the relationship between a drop in voltage and the amount 
of instantaneous electric current in the form of a table, 
to thereby correct the amount of instantaneous electric 
current for each event of the cell. 

Preferably, the modeling step includes 
a step of calculating a drop in voltage from the 
amount of electric current determined according to the 
type of cell, the resistance of a power line, and the 
capacitance of an on-chip capacitor; and 

a correction step of characterizing , for each cell, 
the relationship between a drop in voltage and the amount 
of instantaneous electric current in the form of a table, 
to thereby correct the amount of instantaneous electric 
current for each event of the cell . 

Preferably, the correction step includes a step of 
iterating several times calculation of a drop in voltage 
and correction of a current waveform. 

Preferably, the electromagnetic interference 
analysis method as defined in any one of claims 15 through 
23, wherein the calculation step includes a step of 
modeling the amount of instantaneous electric current 
while separating the same into a short circuit electric 
current component and a charge current component. 
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Preferably, the modeling step includes 
a triangular waveform modeling step of modeling the 
instantaneous current as a triangular waveform whose width 
is calculated for each event information in consideration 
5 of slew information (i.e., an output slew) for an output 
terminal of a cell for each event information such that 
the area of the triangular waveform becomes equal to the 
amount of electric current of each event, the height of 

i£l the triangular waveform being calculated on the basis of 

01 

U 10 the width, and 

fy 

U1 the FFT processing step includes a step of subjecting 

inn 

fi=P to FFT processing information concerning a change in 

current, the information being calculated in the 
™j triangular waveform modeling step. 

^ 15 Preferably, the modeling step includes 

a triangular height calculation step of calculating 
the height of a triangular waveform such that the area of 
the triangular waveform becomes equal to the amount of 
electric current of each event, by means of multiplying 
| 20 the amount of instantaneous electric current by a 

coefficient corresponding to the state of an event of a 
cell, in consideration of whether the event of the cell 
is in a rising state or a falling state. 

Preferably, the modeling step includes 
25 a step of calculating the amount of instantaneous 

26 




electric current for each event information in the case 
of a composite cell; and 

a triangular waveform modeling step of modeling the 
amount of instantaneous electric current as a plurality 
5 of triangular waveforms which are equal in number to the 
stages provided in the composite cell, through use of a 
table representing the relationship between the width and 
height of a triangular waveform; and 
-jg the FFT processing step includes a step of subjecting 

= 10 to FFT processing information concerning a change in 
U1 current, the information being calculated in the 

FU triangular waveform modeling step. 

Further, the present invention provides an 
J 558 ! electromagnetic interference analysis system for 

5f 15 analyzing the amount of electromagnetic interference 

arising in an LSI by means of performing a logic simulation, 
the system comprising: 

a logic simulator; 

computation means which is connected to the logic 
20 simulator and calculates the amount of instantaneous 

electric current from event information, the information 
being produced when a change arises in a signal and 
including the instance name of each cell in which the 
change has arisen, the name of the signal, a time at which 
25 the change has arisen, and transition information; 
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modeling means for modeling the instantaneous 
electric current according to a predetermined rule; and 
fast Fourier (FFT) conversion means for subjecting 
to fast Fourier processing the information concerning a 
5 change in electric current, the information being 

calculated through a modeling step , thereby analyzing the 
amount of electromagnetic interference arising in an LSI 
on the basis of a signal output from the FFT conversion 
means . 

10 

Brief Description of the Drawings 
FIG. 1 is a block diagram showing the overall flow 
of processing pertaining to an EMI analysis method 
according to a first embodiment of the present invention; 
15 FIG. 2 is a block diagram showing a current waveform 

calculation section according to the first embodiment; 

FIG. 3 is a graph showing a calculation result of 
an instantaneous current; 

FIG. 4 is a flowchart showing processing for 
20 averaging a current waveform; 

FIG. 5 is a graph showing a calculation result of 
a current waveform (according to the first embodiment) ; 

FIG. 6 is a graph showing a result of an FFT 
operation; 

25 FIG. 7 shows a rectangular waveform model (according 
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to a second embodiment of the present invention) ; 

FIG, 8 is a block diagram showing calculation of a 
current waveform (according to the second embodiment) ; 

FIG. 9 is a flowchart showing a first rectangular 
waveform-shaping operation ; 

FIG. 10 is a graph showing a calculation result of 
a current waveform (according to the second embodiment) ; 

FIG, 11 shows a rectangular waveform model 
(according to a third embodiment of the present 
invention) ; 

FIG. 12 is a block diagram showing calculation of 
a current waveform (according to the third embodiment) ; 

FIG. 13 is a flowchart showing a second rectangular 
waveform-shaping operation; 

FIG. 14 is a graph showing a calculation result of 
a current waveform (according to the third embodiment) ; 

FIG. 15 shows a rectangular waveform model 
(according to fourth, fifth, and sixth embodiments of the 
present invention) ; 

FIG. 16 shows an i-h table; 

FIG. 17 is a block diagram showing calculation of 
a current waveform (according to a fourth embodiment of 
the present invention) ; 

FIG. 18 is a flowchart showing a third rectangular 
waveform-shaping operation ; 
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FIG. 19 is a graph showing a calculation result of 
a current waveform (according to the fourth embodiment) ; 
FIG. 20 shows an i-s-h table; 

FIG. 21 is a block diagram showing calculation of 
a current waveform (according to a fifth embodiment of the 
present invention) ; 

FIG. 22 is a flowchart showing a fourth rectangular 
waveform-shaping operation ; 

FIG. 23 is a graph showing a calculation result of 
a current waveform (according to the fifth embodiment) ; 

FIG. 24 shows an i-c-h table; 

FIG. 25 is a block diagram showing calculation of 
a current waveform (according to a sixth embodiment of the 
present invention) ; 

FIG. 26 is a flowchart showing a fifth rectangular 
waveform- shaping operation ; 

FIG. 27 is a graph showing a calculation result of 
a current waveform (according to the sixth embodiment) ; 

FIG. 28 is a flowchart for describing an EMI analysis 
method according to a seventh embodiment of the present 
invention ; 

FIG. 29 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
cal cul a ti on/instantaneous current correction section 
according to the seventh embodiment; 
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FIG. 30 is a schematic representation showing layout 
data which are entered into the line voltage drop 
calculation /ins tantaneous current correction section 
according to the seventh embodiment; 
5 FIG. 31 is a schematic representation showing 

information about segments in a power line, which 
information is produced in the line voltage drop 
calculation /ins tantaneous current correction section 
y3 according to the seventh embodiment; 

m 10 FIG. 32 is a schematic representation showing 

y] information about the resistance of each of segments, 

^ which information is produced in the line voltage drop 

*r* calculation/instantaneous current correction section 

?™I according to the seventh embodiment; 

JS 15 FIG. 33 is a schematic representation showing 

information about an electric current flowing in each of 
segments, which information is produced in the line 
voltage drop calculation/instantaneous current 
correction section according to the seventh embodiment; 
20 FIG. 34 is a schematic representation showing 

information about a drop in the voltage of each of segments , 
which information is produced in the line voltage drop 
calculation/instantaneous current correction section 
according to the seventh embodiment; 
25 FIG. 35 is a schematic representation showing 
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information about a drop in the voltage of each instance, 
which information is produced in the line voltage drop 
calculation /instantaneous current correction section 
according to the seventh embodiment; 
5 FIG. 36 is a flowchart for describing an EMI analysis 

method according to an eighth embodiment of the present 
invention ; 

FIG. 37 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
10 calculation/instantaneous current correction section 
according to the eighth embodiment; 

FIG. 38 is a schematic representation showing 
information about a capacitor provided between segments, 
which information is produced in the line voltage drop 
15 calculation/instantaneous current correction section 
according to the eighth embodiment; 

FIG. 39 is a schematic representation showing 
information about a drop in voltage of each of segments, 
which information is produced in the line voltage drop 
20 calculation/instantaneous current correction section 
according to the eighth embodiment; 

FIG. 40 is a flowchart for describing an EMI analysis 
method according to a ninth embodiment of the present 
invention ; 

25 FIG. 41 a flowchart for describing the flow of 
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processing operations of a line voltage drop 
calculation/instantaneous current correction section 
according to a ninth embodiment of the present invention; 

FIG. 42 shows a triangular waveform model (according 
5 to a tenth embodiment of the present invention) ; 

FIG. 43 is a block diagram showing a current waveform 
calculation section (according to a tenth embodiment of 
the present invention) ; 

FIG. 44 is a flowchart of a triangular waveform- 
10 shaping operation (according to the tenth embodiment) ; 

FIG. 45 is a graph showing a calculation result of 
a current waveform (according to the tenth embodiment) ; 

FIG. 46 is a mul ti -order-function model (according 
to an eleventh embodiment of the present invention); 
15 FIG. 4 7 is a block diagram showing a current waveform 

calculation section (according to a tenth embodiment of 
the present invention) ; 

FIG. 48 is a flowchart of a mul ti -order-function 
waveform- shaping operation ; 
20 FIG. 49 is a graph showing a calculation result of 

a current waveform (according to an eleventh embodiment 
of the present invention) ; 

FIG. 50 is a schematic representation showing a 
Gaussian model (according to twelfth and twenty-fifth 
25 embodiments of the present invention) ; 
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FIG. 51 is an i-h table; 
FIG. 52 is an i-w table; 

FIG. 53 is a block diagram showing a current waveform 
calculation section (according to a twelfth embodiment of 
5 the present invention) ; 

FIG. 54 is a flowchart of a Gaussian function 
waveform-shaping operation; 

FIG. 55 is a graph showing a calculation result of 
a current waveform (according to a twelfth embodiment of 
10 the present invention) ; 

FIG. 56 is a block diagram showing the configuration 
of a system for effecting an EMI analysis method according 
to a thirteenth embodiment of the present invention; 

FIG. 57 is an illustration showing example data 
15 stored in means for storing basic information to be used 
for estimating an electric current according to the 
thi rteen th embodiment ; 

FIG. 58 is a graph showing example data stored in 
total current storage means according to the thirteenth 
20 embodiment; 

FIG. 59 is a graph showing example data stored in 
element current storage means according to the thirteenth 
embodiment ; 

FIG. 60 is a flowchart showing processing operations 
25 of the total current calculation means according to the 
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thirteenth embodiment ; 

FIG. 61 is a block diagram showing a current waveform 
calculation section (according to a fourteenth embodiment 
of the present invention) ; 

FIG. 62 is a flowchart of a sixth rectangular 
waveform-shaping operation ; 

FIG. 63 shows an s-h table; 

FIG. 64 is a block diagram showing a current waveform 
calculation section (according to a fifteenth embodiment 
of the present invention) ; 

FIG. 65 is a flowchart of a seventh rectangular 
waveform-shaping operation ; 

FIG. 66 shows a c-h table; 

FIG. 67 is a block diagram showing a current waveform 
calculation section (according to a sixteenth embodiment 
of the present invention) ; 

FIG. 68 is a flowchart of an eighth rectangular 
waveform-shaping operation; 

FIG. 69 is a graph showing example data stored in 
a pattern-of -signal -change storage means according to the 
seventeen thembodimentof the present invention; 

FIG. 70 is a graph showing an example current 
waveform arising in a real circuit according to the 
seventeenth embodiment ; 

FIG. 71 is a graph showing example data stored in 
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conventional total current storage means according to the 
thi r teen th embodiment ; 

FIG. 72 is a block diagram showing the configuration 
of a system for effecting a gate-level EMI analysis method 
5 according to the seventeenth embodiment and an eighteenth 
embodiment of the present invention; 

FIG. 73 is an illustration showing example data 
stored in circuit connection information storage means 

==3 

yp according to the seventeenth embodiment; 

N 5 10 FIG. 74 is a graph showing example data stored in 

y= total current storage means according to the seventeenth 

iy embodiment; 

H FIG. 75 is a graph showing example data stored in 

~*\ element current storage means according to the seventeenth 

5f 15 embodiment; 

FIG. 76 is a flowchart showing processing operations 
of the total current calculation means according to the 
seventeenth embodiment ; 

FIG. 77 is an illustration showing example data 
20 stored in circuit connection information storage means 
according to the eighteenth embodiment; 

FIG. 78 is a graph showing example data stored in 
a pattern-of -signal -change storage means according to the 
eighteenth embodiment ; 
25 FIG. 7 9 is a graph showing an example current 





waveform arising in a real circuit according to the 
eighteenth embodiment ; 

FIG. 80 is a graph showing example data stored in 
conventional total current storage means according to the 
5 eighteenth embodiment; 

FIG. 81 is a graph showing example data stored in 
element current storage means according to the eighteenth 
embodiment ; 

Q 

yB FIG. 82 is a graph showing example data stored in 

M= 10 total current storage means according to the eighteenth 
embodiment; 

FIG. 83 is a flowchart showing processing operations 
of the total current calculation means according to the 
eighteenth embodiment ; 
15 FIG. 84 is an illustration showing example data 

stored in circuit connection information storage means 
according to a nineteenth embodiment of the present 
invention ; 

FIG. 85 is a graph showing example data stored in 
20 a pattern-of -signal -change storage means according to the 
nineteenth embodiment ; 

FIG. 86 is a graph showing an example current 
waveform arising in a real circuit according to the 
nineteenth embodiment ; 
25 FIG. 87 is a graph showing example data stored in 
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conventional total current storage means according to the 
nineteenth embodiment ; 

FIG. 88 is a graph showing example data stored in 
element current storage means according to the nineteenth 
embodiment ; 

FIG . 89 is a flowchart showing processing operations 
of the total current calculation means according to the 
nineteenth embodiment ; 

FIG. 90 is a graph showing example data stored in 
total current storage means according to the nineteenth 
embodiment; 

FIG. 91 is a flowchart for describing an EMI analysis 
method according to a twentieth embodiment of the present 
invention; 

FIG. 92 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
calculation/ instantaneous current correction section 
according to the twentieth embodiment; 

FIG. 93 is a flowchart for describing an EMI analysis 
method according to a twenty-first embodiment of the 
present invention ; 

FIG. 94 is a flowchart for describing the flow of 
processing operations of a line voltage drop 
calculation/instantaneous current correction section 
according to the twenty-first embodiment; 
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FIG. 95 is a flowchart for describing an EMI analysis 
method according to a twenty-second embodiment of the 
present invention ; 

FIG. 96 shows a triangular waveform model (according 
5 to a twenty-third embodiment of the present invention) ; 

FIG. 97 is a block diagram showing a current waveform 
calculation section (according to the twenty-third) ; 

FIG. 98 is a flowchart of a triangular waveform- 
shaping operation ; 
10 FIG. 99 is a multi-order-function model (according 

to a twenty-fourth embodiment of the present invention) ; 

FIG. 100 is a block diagram showing a current 
waveform calculation section (according to the twenty- 
fourth embodiment) ; 
15 FIG. 101 is a flowchart of a multi-order-function 

waveform-shaping operation ; 

FIG. 102 is a block diagram showing a current 
waveform calculation section (according to a twenty-fifth 
embodiment of the present invention) ; 
20 FIG. 103 is a flowchart of a Gaussian function 

waveform- shaping operation ; 

FIG. 104 is a graph showing an example current 
waveform arising in a real circuit according to a 
twenty-sixth embodiment of the present invention; 
25 FIG. 105 is a graph showing example data stored in 
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element current storage means according to the twenty- 
sixth embodiment ; 

FIG. 106 is a graph showing example data stored in 
conventional total current storage means according to the 
5 twenty-sixth embodiment; 

FIG. 107 is a flowchart showing processing 
operations of the total current calculation means 
according to the twenty- sixth embodiment; 

FIG. 108 is a diagram showing example data concerning 
10 internal connections stored in the circuit connection 
information storage means according to the seventeenth 
embodiment shown in FIG. 73; 

FIG. 109 is a block diagram showing the configuration 
of a system for effecting an EMI analysis method employed 
15 in an example conventional EMI analysis method; 

FIG. 110 is an illustration showing example data 
stored in circuit connection information storage means 
employed in an example conventional EMI analysis method 
and in a twenty-sixth embodiment; 
20 FIG. Ill is a graph showing example data stored in 

a pattern-of -signal -change storage means employed in an 
example conventional EMI analysis method and in the 
twenty- sixth embodiment ; 

FIG. 112 is a graph showing example data stored in 
25 element current storage means employed in an example 
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conventional EMI analysis method; 

FIG, 113 is a graph showing example data stored in 
total current storage means employed in an example 
conventional EMI analysis method; 
5 FIG. 114 is a flowchart showing processing 

operations of the total current calculation means employed 
in an example conventional EMI analysis method; 

FIG. 115 is a diagram of changes in signals flowing 

O 

yg through an internal net of a circuit, the changes being 

• 

H 10 detected according to an example conventional EMI analysis 

3*1! i 

U1 method; 

ffk 

py FIG. 116 is a flowchart for describing an example 

conventional EMI analysis method; 
H FIG. 117 is a triangular waveform model according 

S 15 to the present invention (a twenty-seventh embodiment) ; 

FIG. 118 is a block diagram showing calculation of 
a current waveform (according to the twenty-seventh 
embodiment) ; 

FIG. 119 is a flowchart showing shaping of a 
20 triangular waveform (according to the twenty-seventh 
embodiment) ; 

FIG. 120 shows a triangular waveform model according 
to the present invention (a twenty-eighth embodiment) ; 

FIG. 121 is a block diagram showing calculation of 
25 a current waveform (according to the twenty-eighth 
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embodiment) ; 

FIG. 122 shows a triangular waveform model according 
to the present invention (a twenty-ninth embodiment) ; and 

FIG. 123 is a block diagram showing calculation of 
a current waveform (according to the twenty-ninth 
embodiment) . 

Description of the preferred Embodiments ( Firs t 
Embodiment) 

An electromagnetic interference (EMI) analysis 
method according to a first embodiment of the present 
invention will be described hereinbelow. The EMI 
analysis method is characterized by calculating a current 
waveform. 

The method comprises the steps of: 

extracting event information on a per-cell basis 
from a logic simulator; 

calculating from the event information a 
momentarily-changing current of each cell; 

averaging the instantaneously-changing current by 
means of dividing by a certain duration of time, to thereby 
calculate a more realistic change in current; and 

subjecting the calculation result of a current 
waveform to fast Fourier transformation (FFT) , to thereby 
determine a frequency characteristic of electromagnetic 
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interference (EMI) emitted from a circuit of interest, to 
thereby effect EMI analysis. 

FIG. 1 is a block diagram showing the overall flow 
of processing pertaining to the EMI analysis method of the 
5 first embodiment. In the drawing, there are described a 
netlist 1 representing data pertaining to a circuit to be 
subjected to EMI analysis; a test vector 2 for activating 
the netlist 1; a library 3 having described therein 
y3 information concerning individual cells constituting the 

N* 10 netlist 1; a logic simulator 4 for performing a simulation 
|2 on the basis of the information from the library 3; and 

! ~ a current waveform calculation section 7 for calculating 

rl a current waveform from event information 5, the 

rt information having been calculated on a per-cell basis by 

m 15 means of the logic simulator 4 . An FFT processing section 

9 performs an FFT processing operation on the basis of a 
current waveform calculation result 8 calculated by the 
current waveform calculation section 7, to thereby 
calculate an EMI component. Accordingly, there is 
20 obtained an EMI analysis result 10 of the circuit of 

interest. Thus, the EMI analysis method is characterized 
by use of the current waveform calculation section 7 . The 
logic of a cell, delay information concerning the cell, 
and momentary-current-value information concerning the 
25 cell are described in the library 3 according to the types 
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of individual cells . Every time an event arises , the logic 
simulator 4 outputs the event information 5. The term 
"event" used herein designates occurrence of any change 
in a signal within the logic simulator 4. The event 
5 information 5 contains the instance name of an object of 
EMI analysis which accounts for outputting of the event 
information 5, the name of a signal in which change has 
arisen, the time at which the change has arisen, and 
y3 information concerning transition of the signal. The 

N= 10 expression "instance name" corresponds to a proper name 

j£j assigned to each of the cells constituting the netlist 1. 

iy Reference numeral 6 designates a simulation result output 

from the logic simulator 4. 
~1 As can be seen from a detail block diagram shown in 

55 15 FIG. 2, the current waveform calculation section 7 

comprises an instantaneously-changing current 
calculation section 11 and calculates an 
instantaneously-changing current from information 
concerning an instantaneously-changing current of an 
20 instance calculated by the instantaneously-changing 

current calculation section 11, the instance serving as 
an object of EMI analysis. On the basis of a calculation 
result 12 of an instantaneously-changing current 
(hereinafter referred to as an "instantaneously-changing 
25 current calculation result 12"), the current waveform 



calculation section 7 averages a current waveform by means 
of dividing an instantaneously-changing current by a 
predetermined duration of time. The instantaneously- 
changing current calculation section 11 calculates a 
5 momentary change in the current of each cell from the event 
information 5 concerning each cell , the information being 
output from the logic simulator 4. Calculation of 
momentary change in current is technically analogous to 

,s=t 

i—J 

J3 calculation of power consumption through a conventional 

M 10 logic simulation. Upon receipt of the event information 
tTi 5 output from the logic simulator 4, the 

Iy instantaneously-changing current calculation section 11 

extracts, from the library 3, information concerning an 
'~] instantaneously-changing current of an instance, the 

15 instance being an object of EMI analysis, to thereby 
calculate an instantaneously-changing current. 
Reference numeral 12 designates an instantaneously- 
changing current calculation result determined by the 
instantaneously-changing current calculation section 11 . 
20 FIG. 3 is a plot showing an example instantaneously- 
changing current calculation result 12. In the drawing, 
the vertical axis represents the amount of current "it" 
(As), and the horizontal axis represents time "t." 
Reference symbol T designates a duration of time for 
25 averaging a current waveform by means of dividing an 
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instantaneously-changing current by a certain duration of 
time and corresponds to a predetermined constant. At 
designates a minimum time unit to be used for analysis and 
differs from the minimum unit time of the logic simulator 
4 . Reference symbol "t" denotes a current processing time 
and corresponds to a simulation time of the logic simulator 
4 . Reference symbol "it" denotes the value of an 
instantaneously-changing current. As can be seen from 
FIG. 3, according to the instantaneously-changing current 
calculation result, the waveform of an electric current 
assumes a 5 function. Therefore, an EMI component cannot 
be analyzed through an FFT operation. In order to enable 
analysis of an EMI component through FFT, a current 
waveform averaging section 13 averages the 
instantaneously-changing current calculation result 12 by 
means of dividing by a certain duration of time. 

FIG. 4 shows a flowchart of steps relating to 
averaging of a current waveform. The following 
processing operations are repeated until an averaging 
operation is completed; namely, until "t" exceeds an end 
time of simulation (step 1010) . First, a current buffer 
"lb" for averaging purpose is set to 0 (step 1011) . Next, 
i(t) is added to "lb" until the current time "t" becomes 
"t+T" (steps 1012 and 1013). Here, i(t) denotes an 
instantaneously-changing current at time "t" (step 1014) . 
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Next, I ave is assigned to I (x) until a variable *x" shifts 
from -T to "t, " thereby adding At to the variable "x" (steps 
1015 and 1016) . Here, I (x) denotes a current at time "x." 
The current waveform calculation result 8 can be 
5 determined by means of the flowchart of averaging 

operations, FIG. 5 shows a current waveform calculation 
result 8 obtained at time T=5 provided in the 
instantaneously-changing current calculation result 8 

■Q 

J3 shown in FIG. 3. The FFT processing section 9 subjects 

yi! 

H 10 the thus-calculated waveform to an FFT processing 

Jjj operation, thereby determining the frequency 

lu characteristic of an EMI component. In this way, there 

~J can be obtained a result 10 of EMI analysis (hereinafter 

~J referred to as an "EMI analysis result 10") . FIG. 6 shows 

i5 15 a sample of the EMI analysis result 10. From the EMI 

analysis result 10, the frequency characteristic of an EMI 
component arising in the circuit of interest can be 
determined, and the designer of the circuit can study a 
method of addressing EMI . 
20 The EMI analysis method of the present embodiment 

models a current waveform by means of an averaging 
operation and subjecting the thus-modeled waveform to an 
FFT operation, to thereby enable highly-precise analysis 
of EMI. Further, the EMI analysis method of the present 
25 embodiment can analyze an EMI component of a circuit within 



47 



a shorter period of time than can a transistor-level EMI 
analysis method. 
(Second Embodiment) 

An electromagnetic interference (EMI) analysis 
method according to a second embodiment of the present 
invention will be described hereinbelow. The EMI 
analysis method comprises the steps of: 

calculating a momentarily-changing current of each 
cell from the event information 5 output from the logic 
simulator 4 ; 

modeling the thus-calculated momentarily-changing 
current a rectangular waveform having a certain duration 
of time, the peak value of the waveform having been 
calculated such that an instantaneously-changing current 
coincides with the area of the rectangular waveform. 
Accordingly, a more realistic change in current is 
calculated, and the calculation result of a current 
waveform is subjected to FFT operation, wherewith a 
frequency characteristic of an EMI component arising in 
a circuit of interest is determined and EMI analysis is 
effected . 

FIG. 7 shows an instantaneously-changing current 
value and a modeled rectangular waveform according to the 
second embodiment. In the drawing, "i" represents an 
instantaneously-changing current value, and "t" denotes 
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a current processing time. Reference symbol *W" denotes 
the width of a modeled rectangular waveform. The width 
of a rectangular waveforms corresponds to a predetermined 
constant. Further, reference symbol x> h" denotes the 
height of a rectangular waveform, which is to be calculated 
from "i" and "W." 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 8 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously- changing current calculation result 
12" ) determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
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instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
5 Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with that of the current wave calculation 
operation of the first embodiment. The EMI analysis 
y3 method according to the present embodiment comprises a 

M 10 first rectangular waveform height calculation section 20 
for calculating the height of a rectangular waveform such 

t ; ; 

iy that the area of a rectangular waveform serving as a model 

H of the waveform of the instantaneously-changing current 

J~ coincides with the instantaneously-changing current 

S 15 value; and a first rectangular waveform-shaping section 

21 for modeling the instantaneously-changing current 
calculation result 12 as a rectangular waveform having a 
certain duration of time, thereby calculating the waveform 
of an electric current. 
20 FIG. 9 shows a flow of rectangular waveform-shaping 

operations including an operation for calculating the 
height of a rectangular waveform. The following 
operation is repeated on a per-event basis until a current 
wave calculation operation is completed (step 1020) . The 
25 area of a rectangular waveform is defined as x h" and 



assumes a value of The height M h" of the rectangular 

waveform is calculated as i/W (step 1021). Symbol "i" 
denotes an instantaneously-changing current value 
concerning a cell of an event, the event being an object 
of EMI analysis. This processing corresponds to 
processing pertaining to the first rectangular waveform 
height calculation section 20. Next, w h" is added to 
I(t+x) until a variable "x" changes from -W/2 to W/2 , and 
At is added to "x" (steps 1022 and 1023). I (t+x) 
corresponds to the value of electric current flowing 
through all the cells at time t+x. This processing 
corresponds to processing pertaining to the first 
rectangular waveform-shaping section 21. 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
calculation operation and a rectangular waveform-shaping 
operation. FIG. 10 shows an example current wave 
calculation result 8 at W=5 in the instantaneously- 
changing current calculation result shown in FIG. 3. The 
FFT processing section 9 subjects the waveform to an FFT 
operation, to thereby determine the frequency 
characteristic of an EMI component. As a result, there 
can be obtained an EMI analysis result 10. 

According to the EMI analysis method of the present 
embodiment, the waveform of a current is modeled as a 
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rectangular waveform, and the thus -modeled rectangular 
waveform is subjected to an FFT operation, thereby 
enabling more accurate analysis of EMI. The EMI analysis 
method according to the present embodiment enables 
5 processing within a more shorter period of time than does 
the transistor-level EMI analysis method, 
(Third Embodiment) 

An EMI analysis method according to a third 
embodiment of the present invention will now be described. 

10 Under the EMI analysis method according to the present 
embodiment, a momentarily-changing current of each cell 
is calculated from event information concerning to each 
cell, the information being output from a logic simulator. 
The momentarily-changing current is modeled as a 

15 rectangular waveform, wherein a ratio of a duration of time 
to a peak value is made constant such that an 
instantaneously-changing current value becomes equal to 
the area of the rectangular waveform. As a result, there 
is calculated a more realistic change in current, and the 

20 result of calculation of a current waveform is subjected 
to an FFT operation, to thereby determine the frequency 
characteristic of an EMI component of a circuit of interest 
and enable EMI analysis. The EMI analysis method will now 
be described. 

25 FIG. 11 shows a momentary current value and a modeled 
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rectangular waveform according to the third embodiment. 
Reference symbol "i" denotes an instantaneously-changing 
current value, and "t" denotes a current processing time. 
Reference symbol XX K" denotes a ratio of a peak value of 
5 a rectangular waveform to a duration of time. In the 
present embodiment, the ratio is expressed as M peak- 
value/width," The peak-value/width ratio is constant, 
and an appropriate value has been determined beforehand 
as the value of the peak-value/width ratio. Reference 

10 symbol "h" denotes the height of a rectangular waveform, 
and the height M h" is determined from the value of "i" and 
the value of "K." Further, reference symbol "w" denotes 
the width of a rectangular waveform, and the width "w" of 
a rectangular waveform is calculated from the value of "h , " 

15 which has been calculated previously, and the value of "i . " 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

20 FIG. 12 is a detailed block diagram showing the 

current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 

25 from the event information 5 concerning each cell, the 
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information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
5 designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
u instantaneously- changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 

10 instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 8 function. 
Therefore, an EMI component cannot be analyzed through an 

15 FFT operation. Thus far, processing has been identical 
with the current wave calculation operation of the first 
embodiment. 

Reference numeral 30 designates a second rectangular 
waveform height calculation section for calculating the 

20 height of a rectangular waveform such that a ratio of a 
peak value of a rectangular waveform to be modeled to a 
duration of time is constant and such that the area of a 
rectangular waveform becomes equal to an 
instantaneously-changing current value. Reference 

25 numeral 31 designates a second rectangular waveform width 

54 




calculation section for calculating the width of a 
rectangular waveform such that a ratio of a peak value of 
a rectangular waveform serving as a model of the waveform 
of the instantaneously-changing current to a duration of 
5 time is constant and such that the area of a rectangular 
waveform becomes equal to an instantaneously-changing 
current value. Reference numeral 32 designates a second 
rectangular waveform-shaping section for modeling the 
instantaneously-changing current calculation result 12 as 
10 a rectangular waveform having a constant peak-value/width 
ratio. 

FIG. 13 shows the flowchart of processing pertaining 
to the second rectangular waveform-shaping section, 
including processing pertaining to the second rectangular 

15 waveform height calculation section and processing 
pertaining to the second rectangular waveform width 
calculation section. The following processing is 
iterated on a per-event basis until calculation of a 
current waveform is completed (step 1030) . The height of 

20 a rectangular waveform i s expressed by Eq . 1 provided below , 
from the area of a rectangular waveform and = h/w." 
h = iJ(ixK) ... (1) 

The height "h" of a rectangular waveform is calculated 
through use of Eq.l (step 1031). Reference symbol "i" 
25 denotes the value of an instantaneously-changing current 
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of a cell of an event, the event being an object of EMI 
analysis. Reference symbol "K" denotes a ratio of a peak 
value of a rectangular waveform to the width of the same 
and assumes a predetermined constant. This processing 
corresponds to processing pertaining to the second 
rectangular waveform height calculation section 30. The 
area of a rectangular waveform is defined by "w x h . " The 
width "w" of a rectangular waveform is calculated on the 
basis of "w = i/h" (1032) . Reference symbol w h" denotes 
the height xx h" of a rectangular waveform calculated 
through processing pertaining to step 1031. This 
processing corresponds to processing pertaining to the 
second rectangular waveform width calculation section 31 . 
Next, *h" is added to I (t+x) until a variable "x" changes 
from -W/2 to W/2 , and At is added to "x" (steps 1033 and 
1034). I (t+x) corresponds to the value of electric 
current flowing through all the cells at time t+x. This 
processing corresponds to processing pertaining to the 
second rectangular waveform-shaping section 32. 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
calculation operation, a rectangular waveform width 
calculation operation, and a rectangular waveform-shaping 
operation. FIG. 14 shows an example current wave 
calculation result 8 in the instantaneously-changing 
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current calculation result shown in FIG. 3. The FFT 
processing section 9 subjects the waveform to an FFT 
operation, to thereby determine the frequency 
characteristic of an EMI component. As a result, there 
5 can be obtained an EMI analysis result 10. 

According to the EMI analysis method of the present 
embodiment, the waveform of a current is modeled as a 
rectangular waveform having a constant peak-value/width 
ratio, and the thus-modeled rectangular waveform is 

10 subjected to an FFT operation, thereby enabling more 

accurate analysis of EMI . In an LSI circuit involving flow 
of a large amount of electric current, flow of an electric 
current consumes a long time. The EMI analysis method 
according to the present embodiment enables realization 

15 of a more realistic model and hence enables more accurate 
processing. Further, the EMI analysis method enables EMI 
analysis within a much shorter period of time than the 
transistor-level EMI analysis method. 
(Fourth Embodiment) 

20 A fourth embodiment of the present invention will 

now be described. The present embodiment describes a 
method of analyzing EMI by the steps of : 

calculating an instantaneously-changing current of 
each current from event information concerning each cell, 

25 the information being output from a logic simulator; 
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modeling the instantaneously-changing current as a 
rectangular waveform whose area becomes equal to the 
instantaneously-changing current, through use of a global 
table representing the relationship between an 
5 instantaneously-changing current and the height of a 
rectangular waveform, to thereby calculate a more 
realistic change in current; and 

subjecting a calculation result of the current 
waveform to an FFT operation, to thereby determine the 

10 frequency characteristic of an EMI component of a circuit 
of interest. The global table signifies a table which is 
generally prepared without regard to the type of a cell. 

FIG. 15 shows an instantaneously-changing current 
value and a modeled rectangular waveform according to the 

15 present embodiment. Reference symbol xx i" denotes an 

instantaneously-changing current; "t" denotes a current 
processing time; and "h" denotes the height of a 
rectangular waveform serving as a model of the waveform 
of the instantaneously-changing current. The height of 

20 a rectangular waveform is calculated, through use of a 
global table representing the relationship between an 
instantaneously-changing current value and the height of 
a rectangular waveform. Reference symbol "w" denotes the 
width of a rectangular waveform, and the width w w" of a 

25 rectangular waveform is calculated from the value of x> h," 
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which has been calculated previously, and the value of "i . " 

FIG- 16 is a graph showing an example global i-h table 
representing the relationship between an 

instantaneously-changing current value and the height of 
5 a rectangular waveform. The table shows the value of u h" 
at each of certain points The table indicates that 

"h" assumes a value of 0 at i=0, a value of 5 at i=50, and 
a value of 20 at i=100. In reality, the value of "i" 
changes continuously, and hence the table is expressed as 

10 a graph plotting points provided in the table. The value 
of "h" [i.e., h(i)] at a certain value of "i" is expressed 
by Eq.2 provided below, through linear interpolation of 
the values of two points closest to "i" provided in the 
i-h table: one point has lower values (il, hi), and the 

15 other point has higher values (i2, h2 ) . 

h= h2-h\ (/ /1) + M (2) 
/2 - il 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
20 embodiment (see FIG. 1) . 

FIG. 17 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
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calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

In the present embodiment, an i-h table 40 has 
described therein the previously-described global table 
representing the relationship between an 

instantaneously-changing current value and the height of 
a rectangular waveform (see FIG. 16) . A third rectangular 
waveform height calculation section 41 calculates a peak 
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value of* a rectangular waveform serving as a model of an 
instantaneously-changing current, from an 
instantaneously-changing current value and the i-h table 

40. A third rectangular waveform width calculation 

5 section 42 calculates the width of a rectangular waveform 
from a peak value of a rectangular waveform calculated by 
the third rectangular waveform height calculation section 

41, and the instantaneously-changing current value. A 
third rectangular waveform-shaping section 43 models the 

10 instantaneously-changing current calculation result 12 as 
a rectangular waveform, the height of the rectangular 
waveform corresponding to the peak value of the 
rectangular waveform calculated by the third rectangular 
waveform height calculation section 41, and the width of 

15 the rectangular waveform corresponding to the width 
calculated by the third rectangular waveform width 
calculation section 42. 

FIG. 18 is a flowchart showing processing pertaining 
to the third rectangular waveform-shaping section 43, 

20 including processing pertaining to the rectangular 

waveform height calculation section 4 1 and processing 
pertaining to the rectangular waveform width calculation 
section 42. First, the i-h table 40 is read (step 1040) . 
Next, the subsequent processing operations are iterated 

25 on a per-event basis until the current waveform 
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calculation operation is completed (step 1041) . Two 
points closest to "i , " (il , hi) and (i2 , h2) , are extracted 
from the i-h table 40, Reference symbol x> i" is the value 
of an instantaneously-changing current of a cell of an 
5 event, the event being an object of EMI analysis. In the 
table shown in FIG. 16, the two points closest to i=70 are 
(il, hi) = (50, 5) and (i2, h2) = (100, 20) (step 1042). 

The height M h" of a rectangular waveform is 
calculated according to Eq.2 (step 1043) . In connection 

10 with the table shown in FIG. 16, in the case where i=70, 
h(70)=ll. This processing; i.e., processing pertaining 
to steps 1042 and 1043, corresponds to processing 
pertaining to the third rectangular waveform height 
calculation section 41. The area of a rectangular 

15 waveform is defined as "w x h . " The width "w" of a 

rectangular waveform is calculated from w=i/h (step 1044) . 
Reference symbol x> i" denotes the value of an 
instantaneously-changing current of a cell of an event, 
the event being an object of EMI analysis. Reference 

20 symbol xx h" denotes the height of a rectangular waveform 
calculated by the third rectangular waveform height 
calculation section 41. In the table shown in FIG. 16, 
when i=70, w=6 . 3 . This processing corresponds to 
processing pertaining to the third rectangular waveform 

25 width calculation section 42 . Next, "h" is added to I (t+x) 
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until variable u x" changes from -W/2 to W/2 , and At is added 
to variable "x" (steps 1045 and 1046). I (t+x) denotes 
electric currents flowing through all the cells at time 
"t+x." This processing corresponds to processing 
5 pertaining to the third rectangular waveform-shaping 
section 43. 

The current waveform calculation result 8 can be 
determined by means of the flow of a rectangular waveform 
height calculation operation, a rectangular waveform 

10 width calculation operation, and a rectangular 

waveform-shaping operation. FIG. 19 shows an example 
current waveform calculation result 8 obtained in a case 
where the table shown in FIG. 16 is used as the i-h table 
40 and where the instantaneously-changing current value 

15 is 70. The rectangular waveform shown in FIG. 19 is 

subjected to an FFT operation by the FFT processing section 
9, to thereby determine the frequency characteristic of 
an EMI component. Thus, there can be obtained an EMI 
analysis result 10. 

20 (Fifth Embodiment) 

An EMI analysis method according to a fifth 
embodiment of the present invention will now be described. 
The present embodiment describes the EMI analysis method 
comprising the steps of: 

25 calculating an instantaneously-changing current of 
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each cell from event information concerning each cell, the 
event information being output from a logic simulator; 

modeling the instantaneously-changing current into 
a rectangular waveform whose area becomes equal to the 
5 instantaneously-changing current, through use of a global 
table representing the relationship among a slew in an 
input waveform, an instantaneously-changing current value, 
and the height of a rectangular waveform, to thereby 

yy calculate a change in current in consideration of the 

y i 

N 5 10 more-realistic slew in an input waveform; and 

m subjecting the result of calculation of the current 

ry waveform to an FFT operation, thereby determining the 

frequency characteristic of an EMI component arising in 
~i a circuit of interest. 

=S 15 The relationship between the instantaneously- 

changing current value and a modeled rectangular waveform 
according to the present embodiment is the same as in the 
case of the rectangular waveform model according to the 
fourth embodiment (shown in FIG. 15) , except that and 

20 "h" are calculated in consideration of the influence of 
slew in an input waveform. FIG. 20 shows an example global 
i-s-h table representing the relationship among an 
instantaneously-changing current value, the slew in an 
input waveform, and the height of a rectangular waveform, 

25 and a graph plotting the i-s-h table. The table shows the 
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value of "h" at a certain point (i, s) , where *s" denotes 
the slew (or slew) of an input waveform. In connection 
with this table, "h" assumes a value of 0 at i=0 and s=0 ; 
a value of 5 at i=50 and s=0 ; a value of 14 at i=100 and 
s=0 ; a value of 0 at i=0 and s=2 ; a value of 6 at i=50 and 
s=2 ; and a value of 20 at i = 100 and s=2 . In reality, the 
value of "i" and the value of "s" change continuously, and 
hence the table is represented by a three-dimensional 
graph plotting points appearing in the table. The value 
of h(s,i) at a certain value of xx s" and a certain value 
of "i" is expressed by Eq.3 provided below, through use 
of four points in combination, the four points comprising 
a combination of a low value (il) and a high value (i2) 
closest to "i" within the i-s-h table and a low value (si) 
and a high value (s2) closest to "s"; namely, h(sl, il) , 
h(sl, i2) , h(s2, il) , and h(s2, i2). 



(3) 



h(sl 9 /1)02 - s) h(s2, i\)(s - si) 



(52-51) 



f'2-f 
i2~i\J 



h(sl 9 /2)(52 - 5) 7^(52, /2)(5 - 5l) 



(52-51) 



(52-51) 



i-il 
/2-/1 



The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 21 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
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present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
5 information being output from the logic simulator 4 . 

Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an ins tantaneously- 

10 changing current (hereinafter referred to as an 

"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 

15 As can be seen from FIG. 3, according to the 

instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 

20 identical with the current wave calculation operation of 
the first embodiment. 

The i-s-h table 50 employed in the present embodiment 
is the above-described global table representing the 
relationship among an instantaneously-changing current 

25 value, the slew in an input waveform, and the height of 
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a rectangular waveform (FIG. 20) . A fourth rectangular 
waveform height calculation section 52 calculates a peak 
value of a rectangular waveform serving as a model of an 
instantaneously-changing current, from an 
5 instantaneously-changing current value of an instance 
serving as an object of EMI analysis, the slew in an input 
waveform, and the i-s-h table 50, all of which are read 
from Slew data 51. The Slew data 51 are provided in the 
J3 netlist 1 and describe information concerning the slew in 

H 10 a waveform input to each instance. Reference numeral 53 

szr. a 
I U 

yj designates a fourth rectangular waveform width 

r - calculation section for calculating the width of a 

rectangular waveform, through use of the peak value 
~J calculated by the fourth rectangular waveform height 

J= 15 calculation section 52 and the instantaneously-changing 

current value such that the area of the rectangular 
waveform becomes equal to the instantaneously-changing 
current value. Reference numeral 54 designates a fourth 
rectangular waveform-shaping section for modeling the 
20 instantaneously-changing current calculation result 12 as 
a rectangular waveform, the height of the rectangular 
waveform corresponding to the peak value calculated by the 
fourth rectangular waveform height calculation section 52 , 
and the width of the rectangular waveform corresponding 
25 to the width calculated by the fourth rectangular waveform 



67 



width calculation section 53. 

FIG. 22 shows the flowchart of processing pertaining 
to the fourth rectangular waveform-shaping section 54, 
including processing pertaining to the rectangular 
5 waveform height calculation section 52 and processing 
pertaining to the rectangular waveform width calculation 
section 53. First, the i-s-h table 50 is read (step 1050) , 
and the Slew data 51 are read (step 1051) . The processing 
provided below is iterated on a per-event basis until 

10 calculation of a current waveform is completed (step 1052) . 
From the i-s-h table 50, there are extracted four points 
which are combinations of a low value closest to "i," a 
high value closest to "i," a low value closest to "s," and 
a high value closest to "s" ; for example , h (si , il) , h (si , 

15 i2),h(s2 / il), and h (s2 , i2) . Here, reference symbol M i" 
denotes the value of an instantaneously-changing current 
of a cell of an event, the event being an object of EMI 
analysis. Further, reference symbol M s" denotes the slew 
in an input waveform pertaining to the cell. The four 

20 points are extracted on the basis of the name of the 
instance to be processed and the Slew data 51 . 

In connection with the table shown in FIG. 20, in 
the case of i=70 and s=l , there are extracted four values; 
25 h(0, 50)=5, h(0, 100)=14, h(2, 50)=6, and h(2, 100)=20 
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(step 1053) . The height "h" of the rectangular waveform 
is calculated by Eq.3 (step 1054) . In connection with the 
table shown in FIG. 20, in the case where i=70 and s = l, 
h(l, 70) =10.1. Processing pertaining to steps 1053 and 
5 1054 correspond to processing pertaining to the fourth 
rectangular waveform height calculation section 52 . 

The area of the rectangular waveform is defined by 
"w x h . " The duration (w) of a rectangular waveform is 

pi 

J3 calculated on the basis of w=i/h (step 1055) . Here, "i" 

On 

H 6 10 denotes the value of an instantaneously-changing current 
y] of a cell of an event, the event being an object of EMI 

= : 

analysis, and "h" denotes the height of a rectangular 
~ waveform calculated by the fourth rectangular waveform 

ri height calculation section 52 . 

2 15 

In connection with the table shown in FIG. 20, in 
the case of i=70 and s = l , w=6 . 9 . This processing 
corresponds to processing pertaining to the fourth 
rectangular waveform width calculation section 53. Next, 

20 *h" is added to I (t+x) until variable ^x" changes from -W/2 
to W/2 , and At is added to "x" (steps 105 6 and 1057) . Here, 
I (t+x) represents total current flowing through all the 
cells at time "t+x. " This processing corresponds to 
processing pertaining to the fourth rectangular 

25 waveform-shaping section 54 . 



The current waveform calculation result 8 can be 
determined by means of the flow of a rectangular waveform 
height calculation operation, a rectangular waveform 
width calculation operation, and a rectangular 
waveform-shaping operation. FIG. 23 shows an example 
current waveform calculation result 8 obtained in a case 
where the i-s-h table 50 corresponds to the table shown 
in FIG. 20, the instantaneously-changing current value is 
70, and the slew in an input waveform assumes a value of 
1. The rectangular waveform is subjected to an FFT 
operation by the FFT processing section 9, to thereby 
determine the frequency characteristic of an EMI component. 
Thus, there can be obtained an EMI analysis result 10. 
(Sixth Embodiment) 

An EMI analysis method according to a sixth 
embodiment of the present invention will now be described. 
The present embodiment describes the EMI analysis method 
comprising the steps of : 

calculating an instantaneously-changing current of 
each cell from event information concerning each cell, the 
event information being output from a logic simulator; 

modeling the instantaneously-changing current as a 
rectangular waveform whose area becomes equal to the 
instantaneously-changing current, through use of a global 
table representing the relationship among an output load 
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capacitance, an instantaneously-changing current value, 
and the height of a rectangular waveform, to thereby 
calculate a change in current in consideration of the 
more-realistic output load capacitance; and 
5 subjecting the result of calculation of the current 

waveform to an FFT operation, thereby determining the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

The relationship between the ins tantaneously- 

10 changing current value and a modeled rectangular waveform 
according to the present embodiment is the same as in the 
case of the rectangular waveform model according to the 
fourth embodiment (shown in FIG. 15) , except that "w" and 
"h" are calculated in consideration of the influence of 

15 slew in an input waveform. FIG. 24 shows an example global 
i-c-h table representing the relationship among an 
instantaneously-changing current value, an output load 
capacitance, and the height of a rectangular waveform, and 
a graph plotting the i-c-h table. The table shows the 

20 value of M h" at a certain point (i, c) , where "c" denotes 
an output load capacitance . In connection with this table , 
"h" assumes a value of 0 at i=0 and c=0 ; a value of 5 at 
i=5 0 and c=0 ; a value of 20 at i = 100 and c=0 ; a value of 
0 at i = 0 and c=2 0; a value of 8 at i=50 and c=2 0; and a 

25 value of 25 at i = 100 and c=20. In reality, the value of 
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x> i" and the value of w c" change continuously, and hence 
the table is represented by a three-dimensional graph 
plotting points appearing in the table. The value of 
h(c,i) at a certain value of AX c" and a certain value of 
5 is expressed by Eq.4 provided below, through use of 

four points in combination, the four points comprising a 
combination of a low value (il) and a high value (i2) 
closest to "i" within the i-c-h table and a low value (cl) 
and a high value (c2) closest to "c"; namely, h(cl, il), 
10 h(cl, i2) , h(c2, il) , and h(c2, i2). 

y- x + h ... (4) 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 

15 embodiment (see FIG. 1) . 

FIG. 25 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 

20 calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
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conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously- changing current calculation result 
5 12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 

10 waveform of an electric current assumes a 5 function. 

Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

15 The i-c-h table 60 employed in the present embodiment 

is the above-described global table representing the 
relationship between an instantaneously-changing current 
value and the height of a rectangular waveform (FIG. 24) . 
Reference numeral 62 designates a fifth rectangular 

20 waveform height calculation section which calculates a 
peak value of a rectangular waveform serving as a model 
of an instantaneously-changing current, from an 
instantaneously-changing current value of an instance 
serving as an object of EMI analysis and an output load 

25 capacitance. Reference numeral 63 designates a fifth 
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rectangular waveform width calculation section for 
calculating the width of a rectangular waveform, through 
use of the peak value calculated by the fifth rectangular 
waveform height calculation section 62 and the 
5 instantaneously-changing current value and such that the 
area of the rectangular waveform becomes equal to the 
instantaneously-changing current value. Reference 
numeral 64 designates a fifth rectangular waveform- 
shaping section for modeling the instantaneously-changing 

10 current calculation result 12 as a rectangular waveform, 
the height of the rectangular waveform corresponding to 
the peak value calculated by the fifth rectangular 
waveform height calculation section 62, and the width of 
the rectangular waveform corresponding to the width 

15 calculated by the fifth rectangular waveform width 
calculation section 63. 

FIG. 26 shows the flowchart of processing pertaining 
to the fifth rectangular waveform-shaping section 64 , 
including processing pertaining to the rectangular 

20 waveform height calculation section 62 and processing 
pertaining to the rectangular waveform width calculation 
section 63. First, the i-c-h table 60 is read (step 1060) , 
and capacitance data 61 are read (step 1061) . The 
processing provided below is iterated on a per-event basis 

25 until calculation of a current waveform is completed (step 

74 



1062) . From the i-c-h table 60, there are extracted four 
points which are combinations of a low value closest to 
"i," a high value closest to "i," a low value closest to 
"c," and a high value closest to "c"; for example, h(cl, 
5 il), h (cl , i2), h (c2 , il ) , and h (c2 , i2) . Here , reference 
symbol "i" denotes the value of an instantaneously- 
changing current of a cell of an event, the event being 
an object of EMI analysis. Further, reference symbol "c" 
yj denotes an output load capacitance pertaining to a cell 

M= 10 of an event, the event being an object of EMI analysis. 

£ : : 
s *cr 

yi The four points are extracted on the basis of the name of 

instance to be processed and the capacitance data 61 . 

In connection with the table shown in FIG. 24 , in 
~! the case where i=70 and c=10, there are extracted four 

15 values; h(0, 50)=5, h(0, 100)=20, h(20, 50)=8, and h(20, 
100) =25 (step 1063) . The height «h" of the rectangular 
waveform is calculated by Eq . 3 (step 10 64 ) . In connection 
with the table shown in FIG. 24, in the case where i=70 
and 0=10, h(10, 70) =12. 9. Processing pertaining to steps 
20 1063 and 1064 corresponds to processing pertaining to the 
fifth rectangular waveform height calculation section 62 . 
The area of the rectangular waveform is defined by "w x 
h." The duration (w) of a rectangular waveform is 
calculated on the basis of w=i/h (step 1065) . Here, "i" 
25 denotes the value of an instantaneously-changing current 
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of a cell of an event, the event being an object of EMI 
analysis, and "h" denotes the height of a rectangular 
waveform calculated by the fifth rectangular waveform 
height calculation section 62. 

In connection with the table shown in FIG. 24, in 
the case where i=70 and c=10, w=5 . 4 . This processing 
corresponds to processing pertaining to the fifth 
rectangular waveform width calculation section 53. Next, 
"h" is added to I (t+x) until variable "x" changes from -W/2 
to W/2 , and At is added to "x" (steps 1066 and 10 67) . Here, 
I (t+x) represents total current flowing through all the 
cells at time xx t+x." This processing corresponds to 
processing pertaining to the fifth rectangular 
waveform-shaping section 64 . 

The current waveform calculation result 8 can be 
determined by means of the flow of a rectangular waveform 
height calculation operation, a rectangular waveform 
width calculation operation, and a rectangular 
waveform-shaping operation. FIG. 27 shows an example 
current waveform calculation result 8 obtained in a case 
where the i-c-h table 60 corresponds to the table shown 
in FIG. 24, the instantaneously-changing current value is 
70, and an output load capacitance assumes a value of 10. 
The rectangular waveform is subjected to an FFT operation 
by the FFT processing section 9, to thereby determine the 
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frequency characteristic of an EMI component. Thus, 
there can be obtained an EMI analysis result 10. 

The EMI analysis method of the present embodiment 
forms a rectangular waveform model from a global table 
5 defining the relationship between an output load 

capacitance, an instantaneously-changing current value, 
and the height of a rectangular waveform, thereby enabling 
preparation of a more realistic model and more accurate 

O 

y3 processing. 

01 

M 1 10 (Seventh Embodiment) 

Ifi An EMI analysis method according to a seventh 

Iy embodiment of the present invention will now be described. 

H FIG. 28 is a block diagram showing the overall flowchart 

^ of the EMI analysis method according to the present 

15 embodiment. The EMI analysis method according to the 
present embodiment is identical with the EMI analysis 
method according to the first embodiment described in 
connection with FIG. 1, except for addition of data A9 and 
A10, a voltage drop calculation/instantaneously-changing 
20 current correction section All, and an 

instantaneously-changing current calculation result A6 
output in step All. Those portions of the EMI analysis 
method according to the present embodiment which are 
identical with those of the EMI analysis method according 
25 to the first embodiment are assigned the same reference 



numerals, and their repeated explanations are omitted. 
A9 designates data concerning the layout of a 
semiconductor integrated circuit which is an object of EMI 
analysis, A10 designates an I-V table file (or an I-V 
5 function) representing the relationship between the 

amount of electric current and the degree of a drop in line 
voltage. All designates a section of calculating a drop 
in line voltage and correcting the result of an 
instantaneously-changing current calculated in step 11 

10 (hereinafter referred to as a "voltage drop 

calculation/instantaneously-changing current correction 
section") . A6 designates the result of an 
instantaneously-changing current being calculated by the 
vol tage drop cal cul ation /ins tan taneously- changing 

15 current correction section All in consideration of a drop 
in line voltage. 

FIG. 29 is a block diagram showing the flow of 
processing of the voltage drop 

calculation/instantaneously-changing current correction 
20 section All according to the seventh embodiment of the 
present invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 
section All will now be described by reference to FIG. 29. 
Since data 12 have already been described in detail in 
25 connection with the first embodiment, the same reference 
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numeral is assigned to the data 12, and repeated 
explanation is omitted. The voltage drop 

calculation/instantaneously-changing current correction 
section All comprises a power ne^t extraction section B2 , 
a section B6 for calculating the resistance of each segment 
(hereinafter referred to as a ^segment resistance 
calculation section B6") , a section B3 for calculating an 



electric current flowing throx^pjh each segment 
(hereinafter referred to as a "sepient current calculation 

§| 

10 section B3") , a section B7 for 4^jjfc c ulating a drop in line 
voltage of each segment (hereinafter referred to as a 

"segment voltage drop calculati^ section B7") , a section 

* I; 

BIO for calculating a drop in eMne voltage for each 



instance (hereinafter referredpp as an "instance voltage 
15 drop calculation section B10")^*and a current correction 
section B12 . 

These steps will now be described. In step B2 , 

V" 

layout data A9 are entered, ancjf the geometry of a power 

- ■ *. 

line is determined in consideration of information 
20 concerning electrical interconnections. On the basis of 
the thus-determined geometry of/ a power line, nodes, 
change in line width, bends/turns, and switching to 
another interconnection are identified, and the power line 
is branched in accordance with identification results, 
25 thereby forming segment information B4 pertaining to 
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respective segments of the power line. The thus- 

constructed database for the segment information B4 

'v. )■ 

includes the length and width of a line in each segment. 
FIG. 30 is a schematic illustration showing the 
5 layout data A9 which is input to be entered in s.tep B2 . 
FIG. 31 is a schematic illustration showing for th^e data 
base the segment information B4 produced in. step B2. 

Segments 1 through 12 formed as a result of branching of 

\ 

an interconnection are represented by Sgl through S!gl2 . 

10 The segment information B4 concerning respective segments 
of the power line, the information being produced in st^p 
B2 , is entered in step B6, and resistance information B8 X \ 
concerning the resistance of a line in each segment is 
produced from information concerning the length and width 

15 of a line in each segment and a predetermined resistance 
parameter assigned to each interconnection layer. 

FIG. 32 is a schematic representation showing a data 
base for the resistance information B8 produced in step 
B6. In this drawing, the resistance of segment Sgl of a 

20 power line (VDD) is indicated by Rl ; the resistance of 
segment Sg2 of the power line is indicated by R2 ; the 
resistance of segment Sg3 of the power line is indicated 
by R3 ; the resistance of segment Sg4 of the power line is 
indicated by R4 ; the resistance of segment Sg5 of the power 

25 line is indicated by R5 ; and the resistance of segment Sg6 
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of the power line is indicated by R6 . In step B3 , there 
are entered information Bl concerning an 
instantaneously-changing current of each instance 
produced on the basis of the instantaneously-changing 
5 current calculation result 12 (hereinafter referred to 
simply as "instantaneously-changing current information 
Bl") and the segment information B4 concerning the 
respective segments of a power line produced in step B2 . 
The electric current of each instance is distributed into 

10 current components which flow through the respective 

segments of the power line. The current components of each 
instance, which components have been distributed into the 
respective segments of the power line, are summed, to 
thereby produce current information B5 for each segment. 

15 FIG. 33 is a schematic representation showing the 

data base for the current information B5 produced in step 
B3 . In the drawing, electric currents flowing through 
instances 1 through 3 are indicated by II through 13, and 
electrical currents flowing through segments Sgl through 

20 Sg6 are indicated by Isl through Is6. In this case, 

Isl=Il+I2+I3, Is2=I2+I3, Is3=I43, Is4=Il, Is5=I2, and 
Is6=I3. In step B7 , there are entered the current 
information B5 prepared each of the segments in step B3 
and the resistance information B8 prepared for each of the 

25 segments in step B6 . Information B9 concerning a drop in 
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line voltage in each segment (hereinafter referred to 
simply as "voltage drop information B9") is prepared 
according to Ohm's low. 

FIG. 34 is a schematic representation showing the 
data base for the voltage drop information B9 prepared in 
step B7 . In the drawing, voltage drops in respective 
segments Sgl through Sg6 are indicated by AVsl through AVs6 . 
In this case, AVsl=RlxIsl; AVs2=R2xIs2; AVs3=R3xIs3; 
AVs4=R4xIs4; AVs5=R5xIs5; and AVs6=R6xIs6. In step BIO, 
there is entered the voltage drop information B9 prepared 
for each of the segments in step B7 . Voltage drops in 
respective segments of a current flow path of each instance 
are summed, to thereby produce voltage drop information 
Bll for each instance. 

FIG. 35 is a schematic representation showing the 
data base for the voltage drop information Bll prepared 
in step BIO. In the drawing, voltage drops arising in 
instances 1 through 3 are indicated by AVI through AV3 . 
In this case, AVl=AVsl+AVs4 ; AV2=AVsl+AVs2+ AVs5 ; and 
AV3=AVsl + AVs2+AVs3+AVs6 . Here, AVI is expressed as 
AVl=Rlx (11+12+13) +R4xll ; AV2 is expressed as 
AV2=Rlx (11+12+13) + R2x (12+13) +R5xI2 ; and AV3 is expressed 
as AV3=Rlx (I 1+12+ 13) + R2x (12+13) +R3xI3+R6xI3 . 

In step B12 , there are entered the voltage drop 
information Bll prepared for each instance in step BIO and 
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an I-V table file A10 (or an I-V function) representing 
the relationship between the amount of electric current 
and the degree of drop in line voltage, wherewith the 
instantaneously-changing current calculation result 12 is 
corrected- As a result, there is produced an 
instantaneously-changing current calculation result A6 
for each instance, the calculation result taking into 
consideration a drop in line voltage which is dependent 
on the resistance of a power line. 

As mentioned above, the EMI analysis method 
according to the seventh embodiment takes into 
consideration a drop in line voltage which is dependent 
on the resistance of a power line, to thereby correct the 
instantaneously-changing current calculation result 12 
obtained in the first embodiment in consideration of a 
physical phenomena, such as fluctuations in a line voltage 
of a real semiconductor integrated circuit. Therefore, 
the EMI analysis method of the present embodiment can 
produce a more accurate model of a current waveform, 
thereby reducing undesired higher harmonic components 
included in the frequency characteristic of a circuit to 
be analyzed, the characteristic being determined by means 
of subjecting the current waveform to FFT processing. 

Although the seventh embodiment has referred to the 
first embodiment, similar results are yielded even when 
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the technical concept of the present embodiment is applied 
to any of the second through sixth embodiments. 
(Eighth Embodiment) 

An EMI analysis method according to an eighth 
5 embodiment of the present invention will now be described. 
FIG. 36 is a block diagram showing the overall flowchart 
of the EMI analysis method according to the present 
embodiment. The EMI analysis method according to the 
present embodiment is identical with the EMI analysis 
10 method according to the seventh embodiment described in 
connection with FIG. 28, except for addition of data CI 
and replacement of the voltage drop 

calculation/instantaneously-changing current correction 
section All by a voltage drop 

15 calculation/instantaneously-changing current correction 
section C2 . Those portions of the EMI analysis method 
according to the present embodiment which are identical 
with those of the EMI analysis method according to the 
first embodiment are assigned the same reference numerals, 

20 and their repeated explanations are omitted. 

CI designates a C-V table file (or a C-V function) 
representing the relationship between decoupling 
capacitance existing between a power and the ground and 
the degree of a drop in line voltage. C2 designates a 

25 section for calculating a drop in line voltage and 
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correcting the result of an instantaneously-changing 
current being calculated in step 11 (hereinafter referred 
to as a "voltage drop calculation/ins tan taneously- 
changing current correction section") . FIG. 37 is a block 
diagram showing the flow of processing of the voltage drop 
calculation/instantaneously-changing current correction 
section C2 according to the eighth embodiment of the 
present invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 
section C2 will now be described by reference to the 
accompanying drawings . Those portions of the EMI 
analysis method according to the present embodiment which 
are identical with those of the EMI analysis method 
according to the seventh embodiment are assigned the same 
reference numerals, and their repeated explanations are 
omitted. Since the data 12 have already been described 
in detail in connection with the first embodiment, the same 
reference numeral is assigned to the data 12, and its 
repeated explanation is omitted. 

The voltage drop calculation/instantaneously- 
changing current correction section C2 is configured such 
that a section Dl for calculating the capacitance of each 
segment in a power- to-ground line (hereinafter referred 
to simply as a "capacitance calculation section Dl") is 
added to the voltage drop calculation/instantaneously- 
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changing current correction section All, and the segment 
voltage drop calculation section B7 is changed to D3. 
Those sections and their corresponding steps will now be 
described. In step Dl, there is entered the segment 
5 information B4 concerning a power line prepared in step 
B2 in the seventh embodiment. Information D2 concerning 
the capacitance of each segment in the power- to-ground 
line is formed from the information concerning the area 
of an overlap existing between the segments in the 
M= 10 power- to-ground line, the information being included in 

rU 

HI the segment information B4 , as well as from a predetermined 

FU capacitance parameter assigned to each interconnection 

H layer . 

!H FIG. 38 is a schematic illustration showing the data 

15 base for the capacitance information D2 prepared in step 
Dl. In the drawing, a capacitor provided between segments 
Sgl and Sg7 is indicated by CI; a capacitor provided 
between segments Sg2 and Sg8 is indicated by C2 ; and a 
capacitor provided between segments Sg3 and Sg9 is 
20 indicated by C3 . In step D3 , there are entered the current 
information B5 prepared for each segment in step B3 in the 
seventh embodiment, the resistance information B8 
prepared for each segment in step B6 in the seventh 
embodiment, the C-V table file CI (or C-V function) 
25 representing the relationship between the decoupling 
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capacitance existing between the power and the ground and 
the degree of drop in line voltage, and the capacitance 
information D2 prepared in step Dl with regard to the 
capacitance of each of the segments in the power- to-ground 
5 line. Those information pieces are corrected according 
to Ohm's law and the C-V table, to thereby produce the 
voltage drop information B9 for each segment so as to 
reflect decoupling capacitance as well as the resistance 

^0 of a power line. 

01 

M 10 FIG. 39 is a schematic illustration showing the 

voltage drop information B9 prepared in step D3 . In the 

08 

iy drawing, a voltage drop arising in segment Sgl is indicated 

rj by AVsll; a voltage drop arising in segment Sg2 is 

ri indicated by AVs21; a voltage drop arising in segment Sg3 

S5 15 is indicated by AVs31 ; a voltage drop arising in segment 

Sg4 is indicated by AVs41; a voltage drop arising in 
segment Sg5 is indicated by AVsSl; and a voltage drop 
arising in segment Sg is indicated by AVs61 . In this case, 
there are obtained AVsl l=f (CI , AVsl) , AVs21=f (C2, AVs2) , 
20 AVs31=f(C3, AVs3) , AVs41=AVs4, AVs51=AVs5 , and 
AVs61=AVs6 . 

As has been described, the EMI analysis method 
according to the eighth embodiment newly adopts correction 
of a current value in response to fluctuations in a line 
25 voltage, which would be induced by decoupling capacitance 
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existing between the power and the ground, in addition to 
correction of the current value in consideration of" a drop 
in the line voltage, which would be induced by the 
resistance of a power line. The instantaneously-changing 
5 current result 12 obtained for each instance in the first 
embodiment can be corrected in consideration of physical 
phenomena, such as fluctuations in a line voltage of a real 
semiconductor integrated circuit, thereby enabling 
realization of more accurate modeling of a current 

10 waveform. Therefore, the EMI analysis method of the 

present embodiment can produce a more accurate model of 
current waveform, thereby reducing undesired higher 
harmonic components included in the frequency 
characteristic of a circuit to be analyzed, the 

15 characteristic being determined by means of subjecting the 
current waveform to FFT processing. 

Although the eighth embodiment has referred to the 
first embodiment, similar results are yielded even when 
the technical concept of the present embodiment is applied 

20 to any of the second through sixth embodiments . 
(Ninth Embodiment) 

An EMI analysis method according to a ninth 
embodiment of the present invention will now be described. 
FIG. 40 is a block diagram showing the overall flowchart 

25 of the EMI analysis method according to the present 
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embodiment. The EMI analysis method according to the 
present embodiment is identical with the EMI analysis 
method according to the seventh embodiment described in 
connection with FIG. 28, except that the voltage drop 
5 calculation/instantaneously-changing current correction 
section All is replaced by a voltage drop 
calculation/instantaneously-changing current correction 
section El. Those portions of the EMI analysis method 
according to the present embodiment which are identical 

10 with those of the EMI analysis method according to the 
first embodiment are assigned the same reference numerals , 
and their repeated explanations are omitted. El 
designates a section for calculating a drop in line voltage 
and correcting the result of an instantaneously-changing 

15 current calculation performed in step 11 (hereinafter 
referred to as a "voltage drop 

calculation/instantaneously-changing current correction 
section") . 

FIG. 41 is a block diagram showing the flow of 

20 processing of the voltage drop 

calculation/instantaneously-changing current correction 
section El according to the ninth embodiment of the present 
invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 

25 section El will now be described by reference to the 
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accompanying drawings . Those portions of the EMI 
analysis method according to the present embodiment which 
are identical with those of the EMI analysis methods 
according to the seventh and eighth embodiments are 
assigned the same reference numerals, and their repeated 
explanations are omitted. 

Since the data 12 have already been described in 
detail in connection with the first embodiment, the same 
reference numeral is assigned to the data 12, and repeated 
explanation is omitted. The voltage drop 

calculation/instantaneously-changing current correction 
section El is configured by means of removing, from the 
voltage drop calculation /instantaneously- changing 
current correction section C2 employed in the eighth 
embodiment, the segment current calculation section B3 for 
calculating an electric current flowing through each 
segment, the voltage drop calculation section D3 for 
calculating a drop in the voltage of each segment, and the 
instance voltage drop calculation section BIO for 
calculating a drop in the voltage of each instance; and 
newly adding a transient analysis section Fl to the same. 

These sections and their corresponding steps will 
now be described. In step Fl , there are entered the 
instantaneously-changing current information Bl prepared 
for each instance from the instantaneously-changing 

90 




current calculation result 12 in step 11, the resistance 
information B8 prepared for each segment in step B6 in the 
eighth embodiment, and the capacitance information D2 
concerning a capacitor provided between segments prepared 
5 in step Dl in the eight embodiment. Through a simulation 
performed by a transient analysis simulator typified by 
SPICE, there is produced voltage drop information Bll 
concerning a voltage drop of each instance. 
yS As has been described, under the EMI analysis method 

M= 10 according to the ninth embodiment, the resistance of a 
fjf power line and fluctuations in a line voltage, which would 

FU be induced by decoupling capacitance existing between the 

H power line and ground, are determined through use of a 

H transient analysis simulator. Physical phenomena, such 

15 as fluctuations in a line voltage of a real semiconductor 
integrated circuit, can be faithfully reflected n the 
instantaneously-changing current result 12 obtained for 
each instance in the first embodiment. Therefore, the EMI 
analysis method of the present embodiment enables 
20 materialization of highly-accurate EMI analysis as 

described in the background art. Further, there can be 
realized high-speed EMI analysis which can be applied to 
a large-level circuit. 

Although the ninth embodiment has referred to the 
25 first embodiment, similar results are yielded even when 
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the technical idea of the present embodiment is applied 
to any of the second through sixth embodiments. 
(Tenth Embodiment) 

An EMI analysis method according to a tenth 
embodiment of the present invention will now be described. 
In the second through ninth embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 
the information being output from a logic simulator, and 
the thus-calculated instantaneously-changing current is 
modeled as a rectangular waveform, thereby determining 
change in electric current. However, so long as the 
waveform of an instantaneous electric current is modeled 
as a triangular wave, a more realistic change in electric 
current can be determined. By way of example, there will 
now be described an EMI analysis method comprising the 
steps of : 

calculating a peak value such that an 
instantaneously-changing electric current becomes equal 
to the area of a triangular wave; 

modeling an instantaneously-changing electric 
current as a triangular wave having the peak value and a 
certain width, to thereby calculate a more realistic 
change in electric current; and 

subjecting the result of calculation of an electric 
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15 



20 



waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 42 shows an instantaneous electric current and 
a modeled triangular wave according to the present 
embodiment. In the drawing, "i" denotes an 
instantaneously-changing current; and "t" denotes a 
current processing time. Further, "W" denotes the width 
of a triangular wave serving as a model of an 
instantaneously-changing current and corresponds to a 
constant. An appropriate value has been determined as a 
constant beforehand. Reference symbol "h" denotes the 
height of a triangular wave and is calculated on the basis 
of M i" and "W." An electric current at time t+x is 
expressed by Eq.5 provided below. 



h(cj) = 



h(c\J\){c2-c) h(c2,i\)(c-c\) 



(c2 - cl) 



(c2-cl) 



r i2-i 
J2-HJ 



h{c\J2)(c2 - c) h(c2,i2)(c - cl) 



(c2 - cl) 



(c2-c\) 



i-n \ 



... (5) 

Further, an electric current at time t-x is also expressed 
by Eq . 5 . 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 43 is a detailed block diagram showing the 
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current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates. a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 8 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. Reference numeral 100 denotes a 
triangular waveform height calculation section for 
calculating the height of a triangular wave such that the 
area of a triangular wave serving as a model of an 
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instantaneously-changing current becomes equal to an 
instantaneously-changing current. Reference numeral 101 
denotes a triangular-waveform-shaping section for 
modeling the instantaneously-changing current result 12 
5 into a triangular wave having a constant width. 

FIG. 44 is a flowchart of processing pertaining to 
the triangular-waveform-shaping section 101, including 
processing pertaining to the triangular waveform height 
calculation section 100. The following processing is 

10 iterated on a per-event basis until calculation of a 

current waveform is completed (step 1100) . The area of 
a triangular wave is Wxh/2 and assumes a value of an 
instantaneously-changing current "i." The height "h" of 
a triangular wave is calculated from 2xi/W (step 1101). 

15 Reference symbol "i" denotes an ins tantaneously-changing 
current of a cell of an event, the event being an object 
of EMI analysis. This processing corresponds to 
processing pertaining to the triangular waveform height 
calculation section 100. Until variable *x" changes from 

20 0 toW/2, h (c, i) expressed by Eq. 5 is added to I (t+x) and 
I(t-x) . Further, At is added to variable "x" (steps 1102 
and 1103). Here, I (t+x) denotes total electric current 
flowing through all the cells at time t+x, and I(t-x) 
denotes total electric current flowing through all the 

25 cells at time t-x. This processing corresponds to 
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processing pertaining to the triangular-waveform-shaping 
section 101 . 

The current wave calculation result 8 can be obtained 
by means of the flow of a triangular waveform height 
5 calculation operation and a triangular-waveform-shaping 
operation. FIG. 45 shows an example current wave 
calculation result 8 at time W=5 in the 
ins tan taneously- changing current calculation resul t 
shown in FIG. 3. The FFT processing section 9 subjects 

10 the waveform to an FFT operation, to thereby determine the 
frequency characteristic of an EMI component . As a result, 
there can be obtained an EMI analysis result 10. 
(Eleventh Embodiment) 

An EMI analysis method according to an eleventh 

15 embodiment of the present invention will now be described. 
In the second through ninth embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 
the information being output from a logic simulator, and 

20 the thus-calculated instantaneously-changing current is 
modeled as a rectangular waveform, thereby determining 
change in electric current. However, so long as the 
waveform of an instantaneous electric current is modeled 
as a multi -order-function wave, a more realistic change 

25 in electric current can be determined. By way of example, 
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there will now be described an EMI analysis method 
comprising the steps of: 

calculating a peak value such that an 
instantaneously-changing electric current becomes equal 
5 to the area of a triangular wave; 

modeling an instantaneously-changing electric 
current as a multi-order-function having the peak value 
and a given width, to thereby calculate a more realistic 
change in electric current; and 
10 subjecting the result of calculation of an electric 

waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
of interest. 

[0135] 

15 FIG; 46 shows an instantaneous electric current and 

a modeled multi-order-function according to the present 
embodiment. In the drawing, "i" denotes an 
instantaneously-changing current; and "t" denotes a 
current processing time. Further, "W" denotes the width 

20 of a multi-order-function wave serving as a model of an 
instantaneously-changing current and corresponds to a 
constant. An appropriate value has been determined as a 
constant beforehand. Reference symbol M h" denotes the 
height of a multi-order-function wave and is calculated 

25 on the basis of M i" and "W." A multi -order-function wave 
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serving as a model in the present embodiment is expressed 
by Eq.6 provided below. 

y = a(x-W) 2 (x + W) 2 ... (6) . 
The width of a multi-order-function wave model is taken 
as "2 x W. " 

According to Eq.6, the height "h" of a multi- 
order-function wave model is determined when x=0 and is 
expressed by Eq.7 provided below. 

h = axW 4 ... (7) . 

The area of the mul ti -order-function wave model is 
expressed by Eq.8 provided below. 

S = j±axW 5 ... (8) . 

Since the thus-defined area becomes equal to the 
instantaneously-changing current value "i/' "a" is 
expressed though use of M W" and "1" as Eq. 9 provided below. 

a = — ^— (9) . 

\6xJV 5 

As a result of calculation of "a" through use of Eq . 9 , 
the height w h" of a mul ti -order-function wave serving as 
a model is determined. Let W=2 . 5 and i=100, then a=0.96 
is derived from Eq.9. Eq.7 consequently gives h=37.5. 
Further, let W=2 . 5 and i=50, then a=0.4 8 is derived from 
Eq.9. Eq.7 consequently gives h=18.75. 

The overall flow of processing of the EMI analysis 
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method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG, 47 is a detailed block diagram showing the 
5 current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 

O 

y3 from the event information 5 concerning each cell, the 

H> 10 information being output from the logic simulator 4. 
U1 Calculation of momentary change in current is technically 

e 3 S 

fU analogous to calculation of power consumption through a 

JH conventional logic simulation. Reference numeral 12 

~ designates a calculation result of an instantaneously- 

™ 15 changing current (hereinafter referred to as an 

"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. 

FIG. 3 shows an example instantaneously-changing 
20 current calculation result 12. As can be seen from FIG. 
3, according to the instantaneously-changing current 
calculation result, the waveform of an electric current 
assumes a 5 function. Therefore, an EMI component cannot 
be analyzed through an FFT operation. Thus far, the 
25 processing has been identical with the current wave 




calculation operation of the first embodiment . 

Reference numeral 110 denotes a multi-order- 
function waveform height calculation section for 
calculating the height of a multi-order-function wave such 
5 that the area of a multi -order-function wave serving as 
a model of an instantaneously-changing current becomes 
equal to an instantaneously-changing current. Reference 
numeral 111 denotes a multi -order-function waveform- 
shaping section for modeling the instantaneously-changing 
10 current result 12 as a multi-order-function wave having 
a constant width. 
[0142] 

FIG. 48 is a flowchart of processing pertaining to 
the triangular-waveform-shaping section 101, including 

15 processing pertaining to the triangular waveform height 
calculation section 100. The following processing is 
iterated on a per-event basis until calculation of a 
current waveform is completed (step 1110) . A peak- 
value-dependent coefficient "a" of a mul ti -order-function 

20 is calculated from Eq. 9 (step 1111). Reference symbol n i w 
denotes an instantaneously-changing current of a cell of 
an event, the event being an object of EMI analysis, and 
"W" is a constant which is dependent on the width of a wave. 
This processing corresponds to processing pertaining to 

25 the multi-order-function waveform height calculation 
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section 110. 

Until variable "x" changes from -W to W, multi- 
order-function M y" expressed by Eq.6 is added to I(t+x) . 
Further, At is added to variable M x" (steps 1112 and 1113) . 
5 Here, I (t+x) denotes total electric current flowing 
through all the cells at time t+x. This processing 
corresponds to processing pertaining to the multi- 
order-function waveform-shaping section 111. 

The current wave calculation result 8 can be obtained 

10 by means of the flow of a multi-order-function waveform 
height calculation operation and a multi -order-function 
waveform-shaping operation. FIG. 4 9 shows an example 
current wave calculation result 8 at time W=2 . 5 in the 
instantaneously-changing current calculation result 

15 shown in FIG. 3. The FFT processing section 9 subjects 
the waveform to an FFT operation, to thereby determine the 
frequency characteristic of an EMI component . As a result, 
there can be obtained an EMI analysis result 10. 
(Twelfth Embodiment) 

20 An EMI analysis method according to a twelfth 

embodiment of the present invention will now be described. 
In the second through ninth embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 

25 the information being output from a logic simulator, and 
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the thus-calculated instantaneously-changing current is 
modeled as a rectangular waveform, thereby determining 
change in electric current. However, so long as the 
waveform of an instantaneous electric current is modeled 
5 as a Gauss function wave, a more realistic change in 
electric current can be determined. 

By way of example, there will now be described an 
EMI analysis method comprising the steps of: 

calculating the height and width of a Gauss function 
10 wave through use of a global table representing the 
relationship between an instantaneously-changing 
electric current and a Gauss function serving as a model, 
and another table representing the relationship between 
the instantaneously-changing current and the duration of 
15 the Gauss function wave; 

modeling the instantaneously-changing electric 
current as a Gauss function wave, to thereby calculate a 
more realistic change in electric current; and 

subjecting the result of calculation of an electric 
20 waveform to FFT processing, thereby determining the 

frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 50 shows an instantaneous electric current and 
a modeled Gauss function according to the present 
25 embodiment. In the drawing, >x i" denotes an 
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instantaneously-changing current; "t" denotes a current 
processing time; and "h" denotes the height of a Gauss 
function wave serving as a model and calculated from a 
global table representing the relationship between an 
5 instantaneously-changing current and the height of a Gauss 
function serving as a model . 

FIG. 51 is a graph showing an example global i-h table 
representing the relationship between an 
*0 instantaneously-changing current value and the height of 

10 a Gauss function waveform. The table shows the value of 
2 «h" at each of certain points "i." The table indicates 

= y that "h" assumes a value of 0 at i=0, a value of 10 at i=50, 

:~ and a value of 40 at i=100. In reality, the value of "i" 

Jri changes continuously, and hence the table is expressed as 

ii s 

S 15 a graph plotting points provided in the table. The value 

of *h" [i.e., h (i) ] at a certain value of "i" is expressed 
by Eq.2 provided below, through linear interpolation of 
the values of two points closest to "i" provided in the 
i-h table: one point has lower values (il, hi), and the 
20 other point has higher values (i2, h2). Here, "w" 

designates the width of a Gauss function serving as a model , 
and the width "w" of a Gauss function is calculated through 
use of a global table representing the relationship 
between an instantaneously-changing current and the width 
25 of a Gauss function serving as a model. 
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FIG. 52 is a graph showing an example global i-w table 
representing the relationship between an 

instantaneously-changing current value and the width "w" 
of a rectangular waveform. The table shows the value of 
5 "w" at each of certain points "i." The table indicates 
that "w" assumes a value of 0 at 1=0, a value of 1 at 1=20, 
and a value of 1.5 at i=100. In reality, the value of "i" 
changes continuously, and hence the table is expressed as 
a graph plotting points provided in the table . 
10 The value of M w" [i.e. / w(i)] at a certain value of 

"i" is expressed by Eq.10 provided below, through linear 
interpolation of the values of two points closest to "i" 
provided in the i-w table: one point has lower values (il, 
wl) , and the other point has higher values (12., w2 ) . 

15 w(i)= + ... (10) 

/'2-/1 

In the present embodiment, the Gauss function 
serving as a model is expressed by Eq.10. This function 
is called a Gaussian function. Let x=0 , Eq.10 provides 
"h," which represents the height of a Gaussian function 
20 wave. From Eq.10, let x=3xW, there is derived a value of 
0 . OOOlxh (i.e., 0.01% of the height of a Gaussian function 
wave) . This value is taken as a limit on the width of a 
Gaussian function serving as a model, and the width of the 
Gaussian function is taken as "6xw." 
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The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
5 embodiment (see FIG. 1) . 

FIG. 53 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 

10 calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 

15 conventional logic simulation. Reference numeral 12 

designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously -changing current calculation result 
12") determined by the instantaneously-changing current 

20 calculation section 11. 

FIG. 3 shows an example instantaneously-changing 
current calculation result 12. As can be seen from FIG. 
3, according to the instantaneously-changing current 
calculation result, the waveform of an electric current 

25 assumes a 6 function. Therefore, an EMI component cannot 
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be analyzed through an FFT operation. Thus far, the 
processing has been identical with the current wave 
calculation operation of the first embodiment. Reference 
numeral 4 0 designates a global i-h table representing the 
5 relationship between an instantaneously-changing current 
value and the height of a Gaussian function wave (see FIG. 
51) . Reference numeral 120 designates a Gaussian 
function waveform height calculation section for 
calculating, from an instantaneously-changing current 
10 value and the table i-h 40, a peak value of a Gaussian 
function serving as a model . 

Reference numeral 121 designates a global i-w table 
representing the relationship between an 

instantaneously-changing current value and the width of 
15 a Gaussian function wave (see FIG. 52) . Reference numeral 
122 designates a Gaussian function waveform width 
calculation section for calculating the width of a 
Gaussian function wave from the instantaneously-changing 
current value and the i-w table 121. Reference numeral 
20 123 denotes a Gaussian function waveform-shaping section 
for modeling the instantaneously-changing current result 
12 as a Gaussian function wave having a peak value 
calculated by the Gaussian function waveform height 
calculation section 120 and a width calculated by the 
25 Gaussian function waveform width calculation section 122 . 
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FIG. 54 shows the flowchart of processing pertaining 
to the Gaussian function waveform-shaping section 123, 
including processing pertaining to the Gaussian function 
5 waveform height calculation section 120 and processing 
pertaining to the Gaussian function waveform width 
calculation section 122. First, the i-h table 40 is read 
(step 1120) , and the i-w table 121 is read (step 1121) . 
J3 The processing provided below is iterated on a per-event 

H= 10 basis until calculation of a current waveform is completed 

ry 

j£J (step 1122) . There are extracted two points closest to 

iy x> i" provided in the i-h table 40 : one point has lower values 

^ (il, hi) , and the other point has higher values (i2, h2) 

s"""s 

CI (step 1123) . Here, reference symbol "i" denotes the value 

S 15 of an instantaneously-changing current of a cell of an 

event, the event being an object of EMI analysis. 

In connection with the table shown in FIG. 51 , in 
the case of i— 50, there are obtained (il, hl)=(0, 0) and 
(i2, h2) = (50, 5) . In the case of i = 100 , there are obtained 
20 (il, hl) = (50, 5) and (i2, h2) = (100, 20) . FromEq. 2, there 
is calculated the height "h" of the Gaussian function wave 
(step 1124) . In connection with the table shown in FIG. 
51, in the case of i=50, there is obtained h(50)=10. In 
the case of 1=100, there is obtained h(100)=40. This 
25 processing corresponding to processing pertaining to the 
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Gaussian function waveform height calculation section 
120 . 

There are extracted two points closest to xx i" 
provided in the i-w table 121: one point has lower values 
5 (il, wl) , and the other point has higher values (i2, w2) 
(step 1125) . Here, reference symbol "i" denotes the value 
of an instantaneously-changing current of a cell of an 
event, the event being an object of EMI analysis. In 
yj connection with the table shown in FIG. 52, in the case 

M 10 of i=50, there are obtained (il, wl)=(0, 0) and (i2, 

| ==? 

jfj w2) = (20 , 1) . In the case of i=100 , there are obtained (il , 

m wl)={20, 1) and (i2, w2)=(100, 1.5). 

j=j Eq. 10 provides the width v w" of the Gaussian function 

~1 wave (step 1126) . In connection with the table shown in 

^ 15 FIG. 52, in the case of i=50, there is obtained w(50)=1.3. 

In the case of i = 100, there is obtained w(100)=1.5. The 
processing pertaining to steps 1125 and 1126 corresponds 
to processing pertaining to the Gaussian function waveform 
width calculation section 122. 
20 [0158] 

A function expressed by Eq . 11 provided below is added 
to I(t+x) until variable "x" changes from -3xW to 3xW, and 
At is added to u x" (steps 1127 and 1128) . 

y = hxe VwJ ... (11) . 
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Here, I (t+x) represents total current flowing through all 
the cells at time "t+x." This processing corresponds to 
processing pertaining to the Gaussian function 
waveform- shaping section 123. 
5 [0159] 

The current wave calculation result 8 can be obtained 
by means of the flow of a Gaussian function waveform height 
calculation operation, a Gaussian function waveform width 
calculation operation, and a Gaussian function 

10 waveform-shaping operation. FIG. 55 shows an example 

current waveform calculation result 8 as the result of 
calculation of an instantaneously-changing current shown 
in FIG. 3, the calculation result being obtained in a case 
where the table shown in FIG. 51 is used as the i-h table 

15 40 and the table shown in FIG. 52 is used as the i-w table 
121. The current wave calculation result 8 is subjected 
to an FFT operation by the FFT processing section 9, to 
thereby determine the frequency characteristic of an EMI 
component. Thus, there can be obtained an EMI analysis 

20 result 10. 

(Thirteenth Embodiment) 

An electromagnetic interference (EMI) analysis 
method according to a thirteenth embodiment of the present 
invention will now be described. The EMI analysis method 

25 according to the present embodiment is characterized by 
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modeling an instantaneously-changing current while 
information concerning estimation of a current waveform 
(hereinafter referred to simply as "current waveform 
estimation information") is separated into an element 
5 dependent on a short circuit current (also called "short 
circuit-current-dependent element") and another element 
dependent on a charge current (also called "charge- 
current-dependent element") . The expression "short 
circuit current" used herein means an electric current 

10 which flows from the source terminal to the drain terminal 
in each of a p-channel MOS transistor (hereinafter 
abbreviated "p-MOS transistor") and an n-channel MOS 
transistor (hereinafter abbreviated "n-MOS transistor") 
in a case where the p-MOS and n-MOS transistors are 

15 instantaneously and simultaneously brought into an active 
state as a result of a change arising in a signal appearing 
at the gate terminal of a CMOS transistor constituting a 
circuit element. The expression "charge current" means 
an electric current which flows when electric charge is 

20 stored in or discharged from an output load capacitor of 
a circuit element. A short circuit current arises 
substantially simultaneous with occurrence of a change in 
an input signal, and a charge current flows together with 
a change in an output signal at a slight time lag from the 

25 time when change has arisen in an input signal. 
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Accordingly, estimation of a more accurate current 
waveform becomes feasible by means of separate handling 
of the short circuit-current-dependent element of the 
current waveform estimation information and the 
5 charge-current-dependent element . 

FIG. 56 shows the configuration of the EMI analysis 
method according to the present embodiment using a 
technique of averaging a short circuit current and a charge 
current by means of dividing by a discrete width. 

10 A system used with the EMI analysis method shown in 

FIG. 56 has means for analyzing the amount of current 
(hereinafter referred to as w amount-of -current analysis 
means") , the amount-of -current analysis means comprising 
circuit-connection information storage means 7301 for 

15 storing information concerning circuit interconnections; 
pattern-of -signal-change storage means 7302 for storing 
a pattern of signal change; element-current storage means 
7303 for storing the amount of electric current flowing 
through an element ; basic-inf ormation-f or-es timating- 

20 current storage means 7306 for storing basic information 
concerning estimation of an electric current; total- 
amount-of -current storage means 7304 for storing the total 
amount of current flowing through an element; and 
total -amount-of -current calculation means 7305 for 

25 calculating the total amount of electric current. 
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These means constituting the EMI analysis method 
shown in FIG. 56 will now be described. There will now 
be described procedures for calculating the total amount 
of electric current shown in FIG. 58 from circuit 
connection information shown in FIG. 110, a pattern of 
signal change shown in FIG. Ill, and basic information 
concerning estimation of an electric current shown in FIG. 
59. The amount-of -current analysis means comprising 
circuit-connection information storage means 7301 is 
identical with the circuit connection information storage 
means 1001 described in connection with the background art, 
and the pattern-of -signal -change storage means 7302 is 
identical with the signal pattern storage means 1002 
described in connection with the background art . Hence, 
their repeated explanations are omitted for brevity. 

The element-current storage means 7303 stores 
information concerning the amount of electric current 
flowing through an element; that is, information 
concerning the total amount of electric current flowing 
when single change arises in an external terminal of each 
element, while the information is divided into a short 
circuit current component and a charge current component. 
The information concerning the amount of electric current 
flowing through an element, such as that shown in FIG. 57, 
has been stored in the element-current storage means 7303 
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beforehand* Reference numeral 7401 denotes the amount of 
short circuit current which flows when a change arises in 
a signal appearing at a clock input terminal CK of a 
flip-flop FF. Reference numeral 7402 denotes the amount 
5 of charge current which flows when a change arises in a 
signal appearing at a clock input terminal CK of the 
flip-flop FF. Reference numeral 7403 denotes the amount 
of short circuit current which flows when a change arises 
in a signal appearing at an output terminal Q of the 

10 flip-flop FF. Reference numeral 7404 denotes the amount 
of charge current which flows when a change arises in a 
signal appearing at the output terminal Q of the flip- 
flop FF. Reference numeral 7405 denotes the amount of 
short circuit current which flows when a change arises in 

15 a signal appearing at an output terminal Y of a buffer BUF . 
Reference numeral 74 06 denotes the amount of charge 
current which flows when a change arises in a signal 
appearing at the output terminal Y of the buffer BUF. The 
basic-inf ormat ion -for -estimating -cur rent storage means 

20 7306 stores the discrete width of a short circuit current 
and the discrete width of a charge current in the form shown 
in FIG. 59. The basic information concerning estimation 
of an electric current has a discrete width to be used for 
averaging the total amount of electric current which flows 

25 in response to occurrence of a change in a signal, as in 
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the case of the information employed in the first 
embodiment. However, in contrast with the information 
employed in the first embodiment, the basic information 
concerning estimation of an electric current comprises two 
5 pieces of information: that is, information 7601 

concerning the discrete width of a short circuit current, 
and information 7602 concerning the discrete width of a 
charge current. Information 7601 represents the discrete 
y3 width of a short circuit current, and information 7602 

N= 10 represents the discrete width of a charge current. 
j£j The total -amount-of -current calculation means 7305 

' y operates according to a flowchart such as that shown in 

;~j FIG. 60. In contrast with the conventional method, the 

Z\ flowchart shown in FIG. 60 of the EMI analysis method of 

5= 15 the present embodiment is additionally provided with a 

processing operation for averaging the calculated short 
circuit current by means of dividing by the discrete width 
thereof and averaging the calculated charge current by 
means of dividing by the discrete width thereof (step 7706) . 
20 In step 7701, there is read the circuit connection 

information, which is shown in FIG. 110 and is stored in 
the circuit-connection information storage means 7301. 
In step 7702, there is read the pattern of a change in a 
signal , the pattern being stored in the pattern-of- 
25 signal -change storage means 7302. In step 7703, there is 
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read information concerning the amount of electric current 
flowing through an element, the information being stored 
in the element-current storage means 7303 and shown in FIG. 
57 . 

In step 7704, there is read the information 
concerning estimation of a current, the information being 
stored in the basic-inf orma ti on- for-estima ting- cur rent 
storage means 7306 and shown in FIG. 59. In step 7705, 
the pattern of a change in a signal shown in FIG. Ill is 
imparted to the circuit connection information shown in 
FIG. 110, to thereby enable propagation of a signal. In 
a case where a change arises in a signal appearing at an 
external terminal of a circuit element, the element being 
described in the information concerning the amount of 
electric current flowing through an element, a short 
circuit current component and a charge current component 
included in the electric current flowing through the 
external terminal are added, at respective points in time, 
to the short circuit current component and the charge 
current component included in the information concerning 
the total amount of electric current (simply called as 
"total current information"), respectively. 

[0166] 

In step 7706 , the result of an electric current being 
averaged by means of dividing by the discrete width stored 

115 



in the information concerning estimation of a current 
shown in FIG. 59 is added to the short circuit current 
component included in the total current information. The 
result of an electric current being averaged by means of 
dividing by the discrete width stored in the information 
concerning estimation of a current shown in FIG. 59 is 
added to the charge current component included in the total 
current information. In step 7707, after propagation of 
a signal, information concerning the total amount of 
electric current calculated at each points in time shown 
in FIG. 58 is stored in the total -amount-of -cur rent 
storage means 7304. As a result, there is obtained a 
current waveform, such as that shown in FIG. 58, which is 
more realistic than the current waveform obtained by the 
conventional method shown in FIG. 14. 
[0167] 

(Fourteenth Embodiment) 

An EMI analysis method according to a fourteenth 
embodiment of the present invention will now be described. 
The present embodiment describes an EMI analysis method 
comprising the steps of: 

characterizing a peak current of a cell according 
to the type of cell and the amount of instantaneously- 
changing current of each event, to thereby calculate an 
instantaneously-changing current of each cell from event 
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information, the information being output from a logic 
simulator ; 

modeling the instantaneously-changing current as a 
rectangular waveform whose width is calculated such that 
a peak value corresponds to the peak current of the 
rectangular waveform and such that the 

instantaneously-changing current becomes equal to the 
area of the rectangular waveform, to thereby calculate a 
more realistic change in an electric current; and 

subjecting the result of calculation of a current 
waveform to FFT processing, to thereby determine the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

FIG. 15 shows an instantaneously changing current 
and a modeled rectangular waveform according to the 
present embodiment. Here, x \i " denotes an 
instantaneously-changing current, and w t." denotes a 
current processing time. 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 61 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
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instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
"instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

Reference numeral 130 denotes a peak current library. 
A peak current of each cell is characterized according to 
the type of cell, and the thus-characterized value is 
stored in the peak current library 130 . Reference numeral 
131 designates a sixth rectangular waveform height 
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calculation section for extracting, from the peak current 
library 130, a peak value of a rectangular waveform serving 
as a model. Reference numeral 132 designates a sixth 
rectangular waveform width calculation section for 
calculating the width of a rectangular waveform, from the 
peak value extracted by the sixth rectangular waveform 
height calculation section 131 and the 

instantaneously-changing current value, such that the 
instantaneously-changing current becomes equal to the 
area of the rectangular waveform. Reference numeral 133 
denotes a sixth rectangular waveform-shaping section for 
modeling the instantaneously-changing current 
calculation result 12 as a rectangular waveform, the 
height of the rectangular waveform corresponding to the 
peak value extracted by the sixth rectangular waveform 
height calculation section 131, and the width of the 
rectangular waveform corresponding to the width 
calculated by the sixth rectangular waveform width 
calculation section 132. 

FIG. 62 shows the flowchart of processing pertaining 
to the sixth rectangular waveform-shaping section 133, 
including processing pertaining to the sixth rectangular 
waveform height calculation section 131 and processing 
pertaining to the sixth rectangular waveform width 
calculation section 132. First, the peak current library 
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130 is read (step 1130) . The following processing is 
iterated on a per-event basis until calculation of a 
current waveform is completed (step 1130) . Apeak current 
of an instance which is an object of EMI analysis is 
5 extracted from the peak current library 130, and the peak 
current is taken as the height xx h" of a rectangular 
waveform (step 1132) . This processing corresponds to 
processing pertaining to the sixth rectangular waveform 
height calculation section 131. 

10 The area of the rectangular waveform is defined by 

M w x h." The width "w" of a rectangular waveform is 
calculated on the basis of "w = i/h" (step 1133) . 
Reference symbol "i" denotes an instantaneously-changing 
current of a cell of an instance, the instance being an 

15 object of EMI analysis. Reference symbol "h" denotes the 
height "h" of a rectangular waveform calculated through 
processing pertaining to the sixth rectangular waveform 
height calculation section 131 . This processing 
corresponds to processing pertaining to the sixth 

20 rectangular waveform width calculation section 132 . Next , 
"h" is added to I (t+x) until a variable "x" changes from 
-W/2 to W/2 , and At is added to x> x" (steps 1134 and 1135) . 
I (t+x) corresponds to total electric current flowing 
through all the cells at time t+x. This processing 

25 corresponds to processing pertaining to the sixth 
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rectangular waveform-shaping section 133. 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
calculation operation, a rectangular waveform width 
5 calculation operation, and a rectangular waveform-shaping 
operation. The FFT processing section 9 subjects the 
current waveform calculation result 8 to an FFT operation, 
to thereby determine the frequency characteristic of an 
Jj EMI component. As a result, there can be obtained an EMI 

Li I 

Jf; 10 analysis result 10. 

]fj (Fifteenth Embodiment) 

^ An EMI analysis method according to a fifteenth 

s 

H embodiment of the present invention will now be described, 

rj The present embodiment describes an EMI analysis method 

S 15 comprising the steps of : 

characterizing the relationship between a slew in 
a waveform input to each cell and a peak current of the 
cell, according to the type of cell, and the amount of 
instantaneously-changing current for each event, to 
20 thereby calculate an instantaneously-changing current of 
each cell from event information concerning each cell, the 
information being output from a logic simulator; 

modeling the instantaneously-changing current as a 
rectangular waveform whose width is calculated such that 
25 a peak current taking into consideration a slew in the 
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waveform input to the cell corresponds to a peak current 
of the waveform and such that the instantaneously-changing 
current becomes equal to the area of the rectangular 
waveform, to thereby calculate a more realistic change in 
an electric current; and 

subjecting the result of calculation of a current 
waveform to FFT processing, to thereby determine the 
frequency characteristic of an EMI component arising in 
a circuit of interest. 

FIG. 15 shows an instantaneously changing current 
and a modeled rectangular waveform according to the 
present embodiment. Here, "i" denotes an 
instantaneously-changing current, and "t" denotes a 
current processing time. Further, "h" denotes the height 
of a rectangular waveform. The height "h" of a rectangular 
waveform is calculated from a library which characterizes 
the relationship between a slew in the waveform input to 
each cell and the height of a rectangular waveform. 
Reference symbol x> w" denotes the width of a rectangular 
waveform, and the width *w" of a rectangular waveform is 
calculated from the height "h" calculated previously and 
the instantaneously-changing current "i . " The library 
characterizing the relationship between a slew in an input 
waveform and the height of a rectangular waveform is 
prepared according to the type of cell in the form of an 
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s-h table. 

FIG, 63 shows an example s-h table and a graph 
plotting the table. Here, *s" denotes a slew in a waveform 
input to an instance which is an object of EMI analysis. 
5 The table shows the value of "h" at each of certain points 
M s." The table indicates that "h" assumes a value of 0 
at s=0 , a value of 5 at s=l , and a value of 20 at s=2 . In 
reality, the value of "s" changes continuously, and hence 
the table is expressed as a graph plotting points provided 
10 in the table. The value of "h" [i.e., h(s)] at a certain 
value of "s" is expressed by Eq. 12 provided below, through 
linear interpolation of the values of two points closest 
to xx s" provided in the s-h table : one point has lower values 
(si, hi) , and the other point has higher values (s2, h2 ) . 

15 h(s) = h2 ~ hl (5-^l) + M ... (12) 

52 - si 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 
20 FIG. 64 is a detailed block diagram showing the 

current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
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from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
5 conventional logic simulation. Reference numeral 12 

designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
^instantaneously-changing current calculation result 
%Q 12") determined by the instantaneously-changing current 

if 10 calculation section 11. FIG. 3 shows an example 

U] instantaneously-changing current calculation result 12. 

^ y As can be seen from FIG. 3, according to the 

~ instantaneously-changing current calculation result, the 

~1 waveform of an electric current assumes a 5 function. 

52 15 Therefore, an EMI component cannot be analyzed through an 

FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. 

In the present embodiment, an s-h table library 14 0 
20 characterizes the relationship between a slew in a 

waveform input to each cell and the height of a rectangular 
waveform (see FIG. 63) . Reference numeral 51 designates 
Slew data describing a slew in a waveform input to each 
instance. Reference numeral 141 designates a seventh 
25 rectangular waveform height calculation section for 
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extracting a slew in a waveform input to each instance, 
the instance being an object of EMI analysis, from the Slew 
data 51, and calculating a peak value of a rectangular 
waveform serving as a model, from the s-h table library 
5 140 and the thus-extracted slew. Reference numeral 142 
designates a seventh rectangular waveform width 
calculation section for calculating the width of a 
rectangular waveform from the peak value calculated by the 
seventh rectangular waveform height calculation section 

10 141 and from the instantaneously-changing current, such 
that the area of the rectangular waveform becomes equal 
to the instantaneously-changing electric current. A 
seventh rectangular waveform-shaping section 143 models 
the instantaneously-changing current calculation result 

15 12 as a rectangular waveform, the height of the rectangular 
waveform corresponding to the peak value of the 
rectangular waveform calculated by the seventh 
rectangular waveform height calculation section 141, and 
the width of the rectangular waveform corresponding to the 

20 width calculated by the seventh rectangular waveform width 
calculation section 142. 

FIG. 65 is a flowchart showing processing pertaining 
to the seventh rectangular waveform-shaping section 143, 
including processing pertaining to the seventh 

25 rectangular waveform height calculation section 141 and 
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processing pertaining to the seventh rectangular waveform 
width calculation section 142. First, the s-h table 
library 140 is read (step 1140) . The Slew data 51 are read 
(step 1141) . Next, the subsequent processing operations 
are iterated on a per-event basis until the current 
waveform calculation operation is completed (step 1142) . 
Two points closest to "s"; (si, hi) and (s2, h2), are 
extracted from the s-h library table 140. Reference 
symbol "s" denotes a slew in an input waveform, the slew 
being extracted from the Slew data 51 concerning an 
instance which is an object of EMI analysis. 

In connection with the table shown in FIG. 63, in 
the case where s=l . 2 , the two points closest to s=1.2 are 
(si, hi) = (1, 5) and (s2, h2 ) = (2, 20) (step 1143) . The 
height w h" of a rectangular waveform is calculated 
according to Eq. 11 (step 1144) . In the table shown in FIG. 
63, in the case where s=1.2, h(1.2)=8. This processing; 
i.e., processing pertaining to steps 1143 and 1144, 
corresponds to processing pertaining to the seventh 
rectangular waveform height calculation section 141. 

The area of a rectangular waveform is defined as "w 
x h." The width "w" of a rectangular waveform is 
calculated from w=i/h (step 1145) . Reference symbol "i" 
denotes the value of an instantaneously-changing current 
of a cell of an event, the event being an object of EMI 
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analysis. Reference symbol "h" denotes the height of a 
rectangular waveform calculated by the seventh 
rectangular waveform height calculation section 141. 
This processing corresponds to processing pertaining to 
5 the seventh rectangular waveform width calculation 

section 142. Next, xx h" is added to I (t+x) until variable 
"x" changes from -W/2 to W/2, and At is added to variable 
"x" (steps 1146 and 1147) . I (t+x) denotes total electric 
current flowing through all the cells at time 
10 "t+x." This processing corresponds to processing 

pertaining to the seventh rectangular waveform-shaping 
section 143 . 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
15 calculation operation, a rectangular waveform width 

calculation operation, and a rectangular waveform-shaping 
operation. The FFT processing section 9 subjects the 
current waveform calculation result 8 to an FFT operation, 
to thereby determine the frequency characteristic of an 
20 EMI component. As a result, there can be obtained an EMI 
analysis result 10. 
(Sixteenth Embodiment ) 

An EMI analysis method according to a sixteenth 
embodiment of the present invention will now be described. 
25 The present embodiment describes an EMI analysis method 
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comprising the steps of: 

characterizing the relationship between an output 
load capacitance of each cell and a peak current of the 
cell, according to the type of cell, and the amount of 
5 instantaneously-changing current for each event, to 

thereby calculate an instantaneously-changing current of 
each cell from event information concerning each cell, the 
information being output from a logic simulator; 

modeling the instantaneously-changing current as a 

10 rectangular waveform whose width is calculated such that 
a peak current taking into consideration of the output load 
capacitance of the cell corresponds to a peak current of 
the waveform and such that the instantaneously-changing 
current becomes equal to the area of the rectangular 

15 waveform, to thereby calculate a more realistic change in 
an electric current; and 

subjecting the result of calculation of a current 
waveform to FFT processing, to thereby determine the 
frequency characteristic of an EMI component arising in 

20 a circuit of interest. 

FIG. 15 shows an instantaneously changing current 
and a modeled rectangular waveform according to the 
present embodiment. Here, "i" denotes an 
instantaneously-changing current, and "t" denotes a 

25 current processing time. Further, "h" denotes the height 
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of a rectangular waveform. The height "h" of a rectangular 
waveform is calculated from a library which characterizes 
the relationship between the output load capacitance of 
each cell and the height of a rectangular waveform. 
5 Reference symbol "w" denotes the width of a rectangular 
waveform, and the width "w" of a rectangular waveform is 
calculated from the height xx h" calculated previously and 
the instantaneously-changing current "i . " The library 
yg characterizing the relationship between the output load 

H 10 capacitance and the height of a rectangular waveform is 
HI prepared according to the type of cell in the form of a 

^ c-h table . 

H FIG. 66 shows an example c-h table and a graph 

= plotting the table. Here, M s" denotes an output load 

i? 15 capacitance of an instance which is an object of EMI 

analysis. The table shows the value of "h" at each of 
certain points "c." The table indicates that M h" assumes 
a value of 0 at c=0 , a value of 5 at c=10 / and a value of 
20 at c=2 0. In reality, the value of "c" changes 
20 continuously, and hence the table is expressed as a graph 
plotting points provided in the table. The value of "h" 
[i.e., h (c) ] at a certain value of "c" is expressed by Eq. 13 
provided below, through linear interpolation of the values 
of two points closest to "c" provided in the s-h table: 
25 one point has lower values (cl, hi) , and the other point 
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has higher values (c2, h2) . 

h(c)= h2 ~ hX (c- el) + /rl ... (13) 
w c2-cl 

The overall flow of processing pertaining to the EMI 
analysis method of the present embodiment is represented 
5 by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 67 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 

10 instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 

15 analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
* instantaneously -changing current cal culation result 

20 12") determined by the instantaneously-changing current 
calculation section 11. 

FIG. 3 shows an example instantaneously-changing 
current calculation result 12. As can be seen from FIG. 
3, according to the instantaneously-changing current 
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calculation result, the waveform of an electric current 
assumes a 5 function. Therefore, an EMI component cannot 
be analyzed through an FFT operation. Thus far, the 
processing has been identical with the current wave 
5 calculation operation of the first embodiment. Reference 
numeral 150 designates a table library. The c-h table 
library 150 is made up of a c-h table library (see FIG. 
66) characterizing the relationship between the 
kQ previously-described output load capacitance of each cell 

10 and the height of a rectangular waveform. Reference 
■Ul numeral 61 designates capacitance data describing an 

'"y output load capacitance of each instance. Reference 

H numeral 15 J. designates an eighth rectangular waveform 

r? height calculation section for extracting an output load 

z? 15 capacitance of an instance, the instance being an object 

of EMI analysis, from the capacitance data 61 and 
calculating a peak value of a rectangular waveform serving 
as a model, from the c-h table library 150 and the 
thus-extracted output load capacitance. Reference 
20 numeral 152 designates an eighth rectangular waveform 
width calculation section for calculating the width of a 
rectangular waveform from the peak value calculated by the 
eighth rectangular waveform height calculation section 
151 and from the instantaneously-changing current, such 
25 that the area of the rectangular waveform becomes equal 
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to the instantaneously-changing electric current. An 
eighth rectangular waveform-shaping section 153 models 
the instantaneously-changing current calculation result 
12 as a rectangular waveform, the height of the rectangular 
5 waveform corresponding to the peak value of the 

rectangular waveform calculated by the eighth rectangular 
waveform height calculation section 151, and the width of 
the rectangular waveform corresponding to the width 
calculated by the eighth rectangular waveform width 

10 calculation section 152. 

FIG. 68 is a flowchart showing processing pertaining 
to the eighth rectangular waveform-shaping section 153, 
including processing pertaining to the eighth rectangular 
waveform height calculation section 151 and processing 

15 pertaining to the eighth rectangular waveform width 

calculation section 152. First, the c-h table library 150 
is read (step 1150) . The capacitance data 61 are read 
(step 1151) . 

Next, the subsequent processing operations are 
20 iterated on a per-event basis until the current waveform 
calculation operation is completed (step 1152). Two 
points closest to "c"; (cl, hi) and (c2, h2) , are extracted 
from the c-h library table 150. Reference symbol "c" 
denotes a slew in an input waveform, the slew being 
25 extracted from the capacitance data 61 concerning an 
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instance which is an object of EMI analysis. 

In connection with the table shown in FIG. 66, in 
the case where c=12 , the two points closest to c=12 are 
<cl, hi) = (1, 5) and <c2, h2 ) = (20, 20) (step 1153) . The 
height "h" of a rectangular waveform is calculated 
according to Eq. 12 (step 1154) . In the table shown in FIG. 
66 , in the case where 0=12, h (12 ) =8 . This processing ; i . e . , 
processing pertaining to steps 1153 and 1154, corresponds 
to processing pertaining to the eighth rectangular 
waveform height calculation section 151. 

This processing corresponds to processing 
pertaining to the eighth rectangular waveform width 
calculation section 152. Next, "h" is added to I(t+x) 
until variable "x" changes from -W/2 to W/2 , and At is added 
to variable "x" (steps 1156 and 1157). I (t+x) denotes 
total electric current flowing through all the cells at 
time 

"t+x." This processing corresponds to processing 
pertaining to the eighth rectangular waveform-shaping 
section 153 . 

The current wave calculation result 8 can be obtained 
by means of the flow of a rectangular waveform height 
calculation operation, a rectangular waveform width 
calculation operation, and a rectangular waveform-shaping 
operation. The FFT processing section 9 subjects the 
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current waveform calculation result 8 to an FFT operation, 
to thereby determine the frequency characteristic of an 
EMI component. As a result, there can be obtained an EMI 
analysis result 10. 
5 (Seventeenth Embodiment) 

An electromagnetic interference (EMI) analysis 
method according to a seventeenth embodiment of the 
present invention will now be described. Among 
operations related to EMI analysis , the present embodiment 

10 describes particularly extraction of a characteristic of 
EMI and calculation of the amount of an electric current, 
in the case where a circuit element is a composite cell 
formed by combination of a plurality of logic elements. 
In a composite cell, an electric current flows every time 

15 each of logic elements constituting the circuit element 
is activated. A logic element provided in a circuit 
element is activated with a time lag unique to the logic 
element . 

Specifically, when a change arises in an input to 
20 a circuit element, an electric current flows many times 
within the circuit element. In a circuit element 
including logic elements in combination as shown in FIG. 
108, 1N1 is changed with a lag of a gate PI provided in 
a cell (provided that the lag is Ins) from the time changes 
25 have arisen in terminals A and B. Further, Y is changed 
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with a lag of a gate P2 provided in the cell (provided that 
the lag is Ins) from the time a change has arisen in 1N1 . 
In a case where a pattern of signal change shown in FIG. 
69 is imparted to the circuit element (cell) , there arises 
5 a current waveform such as that shown in FIG. 70. 

According to the conventional gate-level EMI 
analysis method, the total amount of electric current 
flowing through each circuit element is considered to be 
dissipated with an involved time lag until a signal appears 

10 at an output pin of the circuit element. For instance, 
in a case where a pattern of signal change shown in FIG. 
69 is imparted to a circuit shown in FIG. 73, the circuit 
employing a circuit element shown in FIG. 108, there is 
output current information such as that shown in FIG. 71. 

15 In contrast, according to the EMI analysis method of the 
present embodiment , an electric current , which would ari se 
when a change occurs in a signal appearing at an input pin 
of a composite cell serving as an object, is registered 
in a library, through use of a peak value of the electric 

20 current and a time lag based on a change in an input signal . 
A logic simulator estimates an electric current, on the 
basis of the assumption that a peak value of an electric 
current includes the previously-described time lag in 
accordance with a change in a signal appearing at an input 

25 pin of a composite cell. 
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FIG. 72 shows the configuration of an EMI analysis 
method using the value of a peak current and the amount 
of electric current. 

A system to be used with the EMI analysis method shown 
5 in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal -change storage 
means 9002; element current storage means 9003 for storing 
the amount of electric current flowing in a circuit 
element; total current storage means 9004 for storing a 

10 total amount of electric current flowing through the 

circuit element; and total current calculation means 9005 
for calculating the total amount of electric current 
flowing through a circuit element. 

Next will be described individual elements 

15 constituting the EMI analysis method shown in FIG. 72 and 
procedures for calculating the total amount of electric 
current shown in FIG. 74 , from the circuit connection 
information shown in FIG. 73, the pattern of signal change 
shown in FIG. 69, and the amount of electric current 

20 flowing through a circuit element shown in FIG. 74. 

The circuit connection information storage means 
9001 stores information concerning the connection of a 
circuit to be analyzed, in the same manner as described 
in connection with the background art. Circuit 

25 connection information, such as that shown in FIG. 73, is 
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stored in advance in the circuit connection information 
storage means 9001. The pattern-of-signal -change 
storage 9002 stores the pattern of signal change to be 
applied to the circuit connection information, as in the 
5 case of a conventional pattern of signal change described 
in connection with the background art. As shown in FIG. 
69, a pattern of change in a signal input to an input port 
CK and a pattern of change in a signal input to an input 
yg port D, such as those shown in FIG. 73, are stored in the 

H 10 pattern-of -signal-change storage 9002 beforehand. 
U1 The element current storage means 9003 stores, for 

each circuit element and in a form such as that shown in 
H FIG. 75, pins 9301 in which electric current arises in an 

™ element after a change has arisen in an external input 

~ 15 terminal of each element, time intervals 9302 (i.e. , the 

time which elapses after a change has arisen in an external 
input terminal), total electric current 9303, and peak 
values 9304. Reference numeral 9305 designates flow of 
an electric current in an ANDOR circuit shown in FIG. 108; 
20 specifically, a total amount of 1 mA including a peak 
current of 0.5 mA flows into the circuit element after 
lapse of 1 ns from when input A has arisen. After lapse 
of 2 ns from when the input A has changed, a total mount 
of 3 mA including a peak value of 1.5 mH flows into the 
25 circuit element. Reference numeral 9306 designates flow 
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of" an electric current in a circuit element shown in FIG. 
108; namely, a total amount of 1 mA including apeak current 
of 0.5 mA flows after lapse of 1 ns from when input B has 
changed. Further, after lapse of 2 ns from when the input 
5 B has changed, an electric current of 3 ma including a peak 
value of 1 . 5 ma flows into the circuit. Reference numeral 
9307 designates flow of an electric current in a circuit 
element shown in FIG. 108; specifically, a total amount 

O 

J3 of 3 mA including a peak value of 1.5 mA flows into the 

M= 10 circuit element after lapse of 2 ns from when input C has 

HI changed . 

rU The total current calculation means 9005 operates 

H according to a flowchart such as that shown in FIG. 76. 

P z According to the flowchart shown in FIG. 7 6 concerning a 

S 15 method of calculating an electric current, processing 

pertaining to step 1504 of the conventional method of 
calculating an electric current is replaced by processing 
(pertaining to step 9404) in which an electric current is 
calculated while a rectangular waveform of the electric 
20 current is assumed on the basis of the total amount of 
electric current and a peak value with a time interval 
provided between information pieces concerning the amount 
of current of an element, the information having been 
stored in the element current storage means in association 
25 with a change in an input signal. In step 9401 is read 
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circuit connection information shown in FIG. 73 stored in 
the circuit connection information storage means 9001. 
In step 9402 is read a pattern of signal change, which is 
shown in FIG. 69 and is stored in the pattern-of- 
5 change-in-signal storage means 9002 . In step 9403 is read 
element current information, which is shown in FIG. 75 and 
is stored in the element current information storage means 
9003. 

yB In step 9404, a pattern of signal change shown in 

M= 10 FIG. 69 is imparted to the circuit connection information 
J£j shown in FIG. 73 in step 9404 , thereby enabling propagation 

iy of a signal . In a case where a change arises in a signal 

™ appearing at an external input terminal of a circuit 

rl element, which element is described in the information 

i 15 concerning the amount of electric current flowing an 

element, a rectangular waveform of an electric current is 
produced, by means of averaging the amount of electric 
current flowing through the external input terminal by 
means of a peak value. The thus-produced rectangular 
20 waveform is added to the information concerning the total 
mount of electric current obtained at a time which lags 
a predetermined period of time from the time a change has 
arisen in the external input terminal . After completion 
of propagation of a signal, in step 9405 the thus- 
25 calculated information concerning the total amount of 
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electric current obtained at each time shown in FIG. 74 
is stored in the total current storage means 9004. As a 
result, there is produced a current waveform, such as that 
shown in FIG. 74/ which is more realistic than that which 
5 is obtained by the conventional current calculation method 
and shown in FIG. 71. 

(Eighteenth Embodiment) 

An electromagnetic interference analysis method 
according to an eighteenth embodiment of the present 
invention will now be described. 

Among operations related to EMI analysis, the 
present embodiment describes in particular extraction of 
a characteristic of EMI and calculation of the amount of 
an electric current, in the case where a circuit element 
is a register (signal holding) element synchronized with 
a clock signal. In a register element such as a flip- 
flop, an electric current flows at every rise or fall of 
a change which arises in a signal appearing at a clock pin. 
For example, a pattern of signal change shown in FIG. 78 
is input to a circuit shown in FIG. 77, an electric current 
shown in FIG. 79 flows through the circuit. As in the case 
of the seventeenth embodiment, an electric current flows 
every time each of the logic elements constituting the 
circuit element is activated. A conventional tool for 
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analyzing power consumption on a gate level represents a 
signal change only when data are input to a flip flop; for 
example, at the time of a rise in a clock signal pin. For 
example, in a case where a pattern of signal change shown 
5 in FIG. 78 is input to the circuit shown in FIG. 77, the 
circuit outputs current information such as that shown in 
FIG. 80. 

In the present embodiment, the timing at which an 
y3 electric current flows through a flip-flop is determined 

M= 10 by means of taking, as a reference, a rise and fall of a 
clock signal appearing at the flip-flop. Therefore, an 
§y electric current is estimated by means of taking, as a 

r"; reference, a rise and fall of the clock signal. FIG. 72 

rt shows the configuration of an EMI analysis method using 

S 15 the value of a peak current and the amount of electric 

current. A system to be used with the EMI analysis method 
shown in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal -change storage 
means 9002; element current storage means 9003 for storing 
20 the amount of electric current flowing in a circuit 

element; total current storage means 9004 for storing a 
total amount of electric current flowing through the 
circuit element; and total current calculation means 9005 
for calculating the total amount of electric current 
25 flowing through a circuit element. 
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[0203] 

Next will be described individual elements 
constituting the EMI analysis method shown in FIG. 72 and 
procedures for calculating the total amount of electric 
5 current shown in FIG. 82, from the circuit connection 
information shown in FIG. 77, the pattern of signal change 
shown in FIG. 78, and the amount of electric current 
flowing through a circuit element shown in FIG. 81. 

Q 

^ The circuit connection information storage means 

gi 

L~ 10 9001 stores information concerning the connection of a 
rL 1 circuit to be analyzed, in the same manner as described 

in connection with the background art. Circuit 
irl connection information, such as that shown in FIG. 77, is 

y stored in advance in the circuit connection information 

2 15 storage means 9001. 

[0205] 

The pattern-of -signal -change storage 9002 stores 
the pattern of signal change to be applied to the circuit 
connection information, as in the case of a conventional - 

20 pattern of signal change described in connection with the 
background art. As shown in FIG. 78, a pattern of change 
in a signal input to an input port CK and a pattern of change 
in a signal input to an input port D, such as those shown 
in FIG. 77, are stored in the pa ttern-of -signal -change 

25 storage 9002 beforehand. 
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The element current storage means 9003 stores, for 
each circuit element and in a form such as that shown in 
FIG. 81, information 9901 concerning the total amount of 
electric current flowing in a circuit element at a rise 
5 and fall of a clock signal and information 9902 about a 
peak, value of the electric current flowing at a rise and 
fall of the clock signal. 

The total current calculation means 9005 operates 
according to a flowchart, such as that shown in FIG. 83. 
10 According to the flowchart relating to a method of 
calculating an electric current shown in FIG. 83, 
processing pertaining to step 1504 of the conventional 
current calculation method is replaced by processing of 
calculating a rectangular waveform of the electric current 
15 from a peak value and a total amount of electric current 
flowing at a rise and fall of a clock signal (step 10104) . 

In step 10101, there is read the circuit connection 
information which is shown in FIG. 77 and stored in the 
circuit connection information storage means 9001. In 
20 step 10102, there is read a pattern of signal change which 
is shown in FIG. 78 and is stored in the pattern-of- 
signal-change storage 9002. 

In step 10103, there is read the element current 
information which is shown in FIG. 114 and stored in the 
25 element current storage means 9003. In step 10104, the 
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pattern of" signal shown in FIG. 78 is imparted to the 
circuit connection information shown in FIG, 77, thereby 
enabling propagation of" a signal . In a case where a change 
arises in a signal appearing at an external input terminal 
5 of a circuit element, which element is described in the 
information concerning the amount of electric current 
flowing an element, a rectangular waveform of an electric 
current is produced by means of averaging an electric 

£0 current corresponding to the electric current flowing 

frt 

M= 10 through the external input terminal at a rise and fall of 

\ y 

W\ the signal by means of a peak value. The thus-produced 

ft? rectangular waveforms are added to the information 

j=f concerning the total mount of electric current obtained 

H| at a time when the change has arisen in the signal appearing 

5= 15 at the external input terminal . 

After completion of propagation of a signal, in step 
10105 the thus-calculated information concerning the 
total amount of electric current obtained at each time 
shown in FIG. 82 is stored in the total current storage 
20 means 9004. As a result, there is produced a current 
waveform, such as that shown in FIG. 82, which is more 
realistic than that obtained by the conventional current 
calculation method and shown in FIG. 80. 
(Nineteenth Embodiment) 
25 An electromagnetic interference analysis method 
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according to a nineteenth embodiment of" the present 
invention will now be described. 

Conventional gate-level power consumption analysis 
means can change total electric current according to the 
5 value of* a signal input to a logic element. More 

specifically, the conventional power consumption analysis 
means enables computation of so-called state-dependent 
power consumption. However, the conventional means 
*e% cannot effect state-dependent power consumption including 

10 secular change in electric current and secular change in 

3 ^ 

HI peak value. 

fy Further, some conventional means determines a 

Q stationary value of another signal in the event of 

C3 occurrence of a change in a certain signal. However, such 

E3 15 means cannot determine a change in another signal . For 

r — ; 

instance, an electric current such as that shown in FIG. 
86 (an electric current may change in amount according to 
whether a change arises in both terminals A and B 
simultaneously or in only terminal A) , on the basis of the 

20 circuit connection information shown in FIG. 84 as well 
as a pattern of signal change shown in FIG. 85. In such 
a case, the conventional means can output only a waveform 
shown in FIG. 87. 

In contrast, the EMI analysis method according to 

25 the present embodiment is characterized by holding secular 
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change in electric current and a peak value in the form 
of a state-dependent signal which takes into consideration 
changes in two or more signals. FIG. 72 shows the, 
configuration of an embodiment of the EMI analysis method 
5 using the peak current and the amount of electric current, 
which are calculated in the fourteenth embodiment. 

A system to be used with the EMI analysis method shown 
in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal -change storage 
Uz 10 9002; element current storage means 9003 for storing the 

m amount of electric current flowing in a circuit element; 

83 

fU total current storage means 9004 for storing a total amount 

Q of electric current flowing through the circuit element; 

■3. E 

H and total current calculation means 9005 for calculating 

■y 15 the total amount of electric current flowing through a 

S ,. Ml j 

circuit element. 

Next will be described individual elements 
constituting the EMI analysis method shown in FIG. 72 and 
procedures for calculating the total amount of electric 
20 current shown in FIG. 90, from the circuit connection 
information shown in FIG. 84, the pattern of signal change 
shown in FIG. 85, and the amount of electric current 
flowing through a circuit element shown in FIG. 88. 

The circuit connection information storage means 
25 9001 stores information concerning the connection of a 
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circuit to be analyzed, in the same manner as described 
in connection with the background art. Circuit 
connection information, such as that shown in FIG. 84, is 
stored in advance in the circuit connection information 
5 storage means 9001 . 

The pattern-of -signal -change storage 9002 stores 
the pattern of signal change to be applied to the circuit 
connection information, as in the case of a conventional 
pattern of signal change described in connection with the 

10 background art. As shown in FIG. 85, a pattern of change 
in a signal input to an input port A and a pattern of change 
in a signal input to an input port B, such as those shown 
in FIG. 84, are stored in the pattern-of -signal -change 
storage 9002 beforehand. 

15 The element current storage means 9003 stores, for 

each circuit element and in a form such as that shown in 
FIG . 88 , a changed pin concerning an electric current which 
flows in a circuit element in the event of a change arising 
in the circuit element when a change occurs in an external 

20 terminal of the circuit element, a change arising in the 
logic status of another external pin under the influence 
of the change in the external terminal, total electric 
current, and information concerning a peak value. 

25 For example, reference numeral 10601 designates a 
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circuit element: that is, a NAND circuit shown in FIG. 84, 
in which a total amount of 3 mA including a peak current 
of 1 mA flows into the circuit element, provided that an 
input signal A assumes a logic value of 1 and the logic 
value of an input signal B changes from 0 to 1 when a change 
arises in an output Y. 

The total current calculation means 9005 operates 
according to a flowchart, such as that shown in FIG. 89. 
According to the flowchart relating to a method of 
calculating an electric current shown in FIG. 89, 
processing pertaining to step 1504 of the conventional 
current calculation method is replaced by processing of 
calculating a rectangular waveform of the electric current 
from total electric current and a peak value at the time 
of occurrence of a change in consideration of the status 
of another terminal (step 10704) . 

In step 10701, there is read the circuit connection 
information which is shown in FIG. 84 and stored in the 
circuit connection information storage means 9001. In 
step 10702, there is read a pattern of signal change which 
is shown in FIG. 85 and is stored in the pattern-of- 
signal-change storage 9002. In step 10703, there is read 
the element current information which is shown in FIG. 88 
and stored in the element current storage means 9003. 
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In step 10704, the pattern of signal shown in FIG. 
85 is imparted to the circuit connection information shown 
in FIG. 84, thereby enabling propagation of a signal. In 
a case where a change arises in a signal appearing at an 
5 external input terminal of a circuit element, which 

element is described in the information concerning the 
amount of electric current flowing in an element, a 
rectangular waveform of an electric current is produced 
by means of averaging the amount of electric current 

10 flowing through the external input terminal by means of 
a peak value in accordance with the status of another 
terminal. The thus -produced rectangular waveform is 
added to the information concerning the total mount of 
electric current obtained at a time when the change has 

15 arisen in the signal appearing at the external input 
terminal . 

After completion of propagation of a signal, in step 
10705 the thus-calculated information concerning the 
total amount of electric current obtained at each time 

20 shown in FIG. 90 is stored in the total current storage 
means 9004. As a result, there is produced a current 
waveform, such as that shown in FIG. 90, which is more 
realistic than that which is obtained by the conventional 
current calculation method and shown in FIG. 87. 

25 (Twentieth Embodiment) 
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FIG. 91 is a block diagram showing the overall 
flowchart of an EMI analysis method according to the 
present embodiment. The EMI analysis method according to 
the present embodiment is identical with the EMI analysis 
5 method according to the seventh embodiment described in 
connection with FIG. 28, except for addition of a cell 
characterization section G2 , replacement of the I-V table 
file A10 by an I-V table file (or an I-V function) G3 , and 
replacement of the voltage drop 
10 calculation/instantaneously-changing current correction 
section All by a line voltage drop 

calculation/instantaneously-changing current correction 
section G4 . 

Reference symbol G2 denotes a cell characterization 
15 section for producing information unique to each cell; G3 
denotes an I-V table file (or an I-V function) which 
represents, for each cell, the relationship between the 
amount of electric current and the degree of drop in line 
voltage; and G4 denotes a line voltage drop 
20 calculation/instantaneously-changing current correction 
section for calculating the degree of drop in line voltage 
and correcting the instantaneously-changing current 
calculated in step 11. 
[0226] 

25 The operation of the cell characterization section 
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G2 will now be described. The cell characterization 
section G2 receives a library 3 and calculates in advance 
the amount of electric current at the time of a change 
arising in a line voltage of each cell, thereby producing 
5 an I-V table file (or an I-V function) for each cell. 

FIG. 92 is a block diagram showing the flow of 
processing of the voltage drop 

calculation/instantaneously-changing current correction 
section G4 according to the twentieth embodiment of the 

10 present invention. The function of the voltage drop 

calculation/instantaneously-changing current correction 
section G4 will now be described by reference to FIG. 92. 

Those portions of the EMI analysis method according 
to the present embodiment which are identical with those 

15 of the EMI analysis method according to the seventh 

embodiment are assigned the same reference numerals, and 
their repeated explanations are omitted. Further, since 
the data 12 have already been described in detail in 
connection with the first embodiment, the same reference 

20 numeral is assigned to the data 12, and repeated 

explanation is omitted. In contrast with the voltage drop 
calculation/instantaneously-changing current correction 
section All employed in the seventh embodiment, the 
voltage drop calculation /instantaneously-changing 

25 current correction section G4 is configured such that the 
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current value correction section B12 is replaced by step 
HI . 



There will now be described step HI. In step HI, 
there are entered the voltage drop information Bll 
5 prepared for each instance in step BIO according to the 
seventh embodiment and the I-V table file (or I-V function) 
which is produced for each cell by the cell 
characterization section G2 and represents the 

D 

yQ relationship between the amount of electric current and 

10 the degree of drop in line voltage. On the basis of a table 
yj file unique to each cell , the instantaneously-changing 

current information Bl , which is produced on the basis of 
H the instantaneously-changing current result 12 produced 

H: in step 11, is corrected. An instantaneously-changing 

5? 15 current calculation result A6 is produced for each 

instance in consideration of a drop in line voltage which 
is dependent on the resistance of a power line. 



As mentioned above, under the EMI analysis method 
20 according to the twentieth embodiment of the present 

invention , the relationship between the amount of electric 
current and the degree of drop in line voltage is 
characterized for each cell beforehand. In the seventh 
embodiment, EMI of a circuit element is analyzed on the 
25 basis of a uniform relationship between the amount of 
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electric . current and the degree of drop in line voltage 
using the I-V table file (or I-V function) . In contrast 
with the EMI analysis method of the seventh embodiment, 
the EMI analysis method of the present embodiment enables 
5 dependence of EMI on the structure (or layout) of a cell 
to be reflected in EMI analysis. Therefore, the EMI 
analysis method of the present embodiment can produce a 
more accurate model of a current waveform, thereby 
reducing undesired higher harmonic components included in 

10 the frequency characteristic of a circuit to be analyzed, 
the characteristic being determined by means of subjecting 
the current waveform to FFT processing. 

Although the twentieth embodiment has described by 
reference to the first embodiment, the same result will 

15 be yielded even when the technical idea of the present 
embodiment is applied to any one of the fourteenth through 
nineteenth embodiments . 
(Twenty-first Embodiment) 

FIG. 93 is a block diagram showing the overall 

20 flowchart of an EMI analysis method according to a 

twenty-first embodiment of the present invention. The 
EMI analysis method according to the present embodiment 
is identical with the EMI analysis method according to the 
eighth embodiment described in connection with FIG. 36, 

25 except for addition of a cell characterization section II, 
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deletion of the C-V table file (or a C-V function) CI, 
replacement of the I-V table file A10 by an I -C-V table 
file (or an I-V function) 12, and replacement of the 
vol tage drop calculation/ instantaneously-changing 
current correction section C2 by a line voltage drop 
calculation/instantaneously-changing current correction 
section 13. 

Those portions of the EMI analysis method according 
to the present embodiment which are identical with those 
of the EMI analysis method according to the first 
embodiment are assigned the same reference numerals, and 
their repeated explanations are omitted. Reference 
symbol II denotes a cell characterization section for 
producing information unique to each cell; 12 denotes an 
I-C-V table file (or an I-C-V function) which represents, 
for each cell, the relationship among the amount of 
electric current, the degree of drop in line voltage, and 
decoupling capacitance existing between a power line and 
ground; and 13 denotes a line voltage drop 
calculation/instantaneously-changing current correction 
section for calculating the degree of drop in line voltage 
and correcting the instantaneously-changing current 
calculation result 12 produced in step 11. 

The operation of the cell characterization section 
II will now be described. The cell characterization 
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section II receives a library 3, changes a line voltage 
of each cell, and calculates in advance the amount of 
electric current at the time of changing of the line 
voltage in consideration of decoupling capacitance 
5 existing between a power and ground, thereby producing an 
I-C-V table file (or an I-C-V function) 12 for each cell. 

FIG. 94 is a block diagram showing the flow of 
processing of the voltage drop 

calculation/instantaneously-changing current correction 

10 section 14 according to the twenty-first embodiment of the 
present invention. The function of the voltage drop 
calculation/instantaneously-changing current correction 
section 13 will now be described by reference to FIG. 94. 
Those portions of the EMI analysis method according to the 

15 present embodiment which are identical with those of the 
EMI analysis method according to the eighth embodiment are 
assigned the same reference numerals, and their repeated 
explanations are omitted. Further, since the data 12 have 
already been described in detail in connection with the 

20 first embodiment, the same reference numeral is assigned 
to the data 12, and repeated explanation is omitted. 

In contrast with the voltage drop 
calculation/instantaneously-changing current correction 
section C2 employed in the eighth embodiment, the voltage 

25 drop calculation /instantaneously-changing current 
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correction section 13 is configured such that the current 
value correction section B12 is replaced by step Jl . 

There will now be described step Jl . In step 
Jl , there are entered the voltage drop information Bll 
5 prepared for each instance in step BIO according to the 
eighth embodiment, the I-C-V table file (or I-C-V 
function) 12 which is produced for each cell by the cell 
characterization section II and represents the 
relationship between the amount of electric current and 

10 the degree of drop in line voltage, and the capacitance 
information D2 concerning a capaci tor provided in a 
segment interconnecting a power and ground, the 
information being prepared in step Dl in the eight 
embodiment. On the basis of a table file unique to each 

15 cell , the instantaneously-changing current information Bl , 
which is produced on the basis of the instantaneously- 
changing current result 12 produced in step 11, is 
corrected. An instantaneously-changing current 
calculation result A6 is produced for each instance, in 

20 consideration of fluctuations in line voltage which are 
dependent on decoupling capacitance as well as the 
resistance of a power line. 

As mentioned above, under the EMI analysis method 
according to the twenty-first embodiment of the present 

25 invention, the relationship among the amount of electric 
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current, the degree of drop in line voltage, and decoupling 
capacitance existing between a power and ground is 
characterized for each cell beforehand. In the eighth 
embodiment, EMI of a circuit element is analyzed on the 
5 basis of a uniform relationship between the amount of 
electric current, the degree of drop in line voltage, and 
decoupling capacitance existing between a power and ground, 
using the I-V table file (or I-V function) A10 and the C-V 
table (or C-V function) CI. In contrast with the EMI 

10 analysis method of the eighth embodiment, the EMI analysis 
method of the present embodiment enables dependence of EMI 
on the structure (or layout) of a cell to be reflected in 
EMI analysis. Therefore, the EMI analysis method of the 
present embodiment can produce a more accurate model of 

15 a current waveform, thereby reducing undesired higher 
harmonic components included in the frequency 
characteristic of a circuit to be analyzed, the 
characteristic being determined by means of subjecting the 
current waveform to FFT processing. 

20 Although the twenty-first embodiment has described 

by reference to the first embodiment, the same result will 
be yielded even when the technical idea of the present 
embodiment is applied to any one of the fourteenth through 
nineteenth embodiments . 

25 (Twenty-second Embodiment) 
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FIG. 95 is a block diagram showing the overall 
flowchart of an EMI analysis method according to a 
twenty-second embodiment of the present invention. The 
EMI analysis method according to the present embodiment 
5 is identical with the EMI analysis method according to the 
seventh embodiment described in connection with FIG. 28, 
except for addition of a compari son-and-de termination 
section Kl . Those portions of the EMI analysis method 
according to the present embodiment which are identical 

m 

10 with those of the EMI analysis method according to the 
Ln first embodiment are assigned the same reference numerals, 

fy and their repeated explanations are omitted. Reference 

s 

O symbol Kl denotes a compari son-and-de termination section 

w for determining whether or not the instantaneously- 

hf 15 changing current calculation result A6 — which takes into 
consideration a drop in line voltage produced by the 
vol tage drop calculation/ins tan taneously- changing 
current correction section All — becomes equal to an input 
to the next processing step. 
20 Next will be described the operation of the 

comparison-and-determination section Kl . The EMI 
analysis method according to the present embodiment is 
characterized by means of iterating the following 
operations; that is, re-entry, into the voltage drop 
25 calculation/instantaneously-changing current correction 
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section All, of the instantaneously-changing current 
calculation result A6 which takes into consideration a 
drop in line voltage produced by the voltage drop 
calculation/instantaneously-changing current correction 
5 section All, thereby producing another 

instantaneously-changing current calculation result A6 
which newly takes into consideration a drop in line 
voltage. 

yg The comparison-and-de termination section Kl 

M= 10 determines, through comparison, whether or not the number 
U1 of times the foregoing operations has been iterated has 

RJ reached a predetermined number of times. If a truth-value 

H TRUE is yielded as a result of comparison, iteration of 

^ the operations is completed, and there is performed 

JS 15 processing pertaining to the next step; that is, the 

operation of the FFT processing section 9, is performed. 

As mentioned above, under the EMI analysis 
method according to the twenty-second embodiment of the 
present invention, the instantaneously-changing electric 
20 current result which has been corrected and calculated in 
consideration of the influence of fluctuations in line 
voltage is fed back, to thereby iterate voltage drop 
calculation/instantaneously-changing current correction 
operations. Accordingly, variation in the amount of 
25 electric current, which would otherwise be caused by a drop 




in line voltage, can be made stable. Therefore, the EMI 
analysis method of the present embodiment can produce a 
more accurate model of a current waveform, thereby 
reducing undesired higher harmonic components included in 
5 the frequency characteristic of a circuit to be analyzed, 
the characteristic being determined by means of subjecting 
the current waveform to FFT processing. 

In the present embodiment, the comparison-and- 
g3 determination section Kl controls iteration of processing 

H 1 10 operations by means of determining whether or not a 

Ul predetermined number of times has been reached. 

f« Alternatively, there may be employed a storage device for 

temporarily storing the instantaneously-changing current 
H calculation result A6 which takes into consideration a 

y 15 drop in line voltage produced by the voltage drop 

calculation/instantaneously-changing current correction 
section All. The instantaneously-changing current 
calculation result A6 stored in the storage device is 
entered again into the voltage drop 
20 calculation/instantaneously-changing current correction 
section All, to thereby newly produce another 
instantaneously-changing current calculation result A6 . 
Comparison is made between the thus-produced another 
instantaneously-changing current calculation result A6 
25 and the instantaneously-changing current calculation 




result A6 stored in the storage device. The foregoing 
operations are iterated until a difference falls within 
a predetermined range of allowance. 

Although the twenty-second embodiment has been 
5 described with reference to the seventh embodiment, the 
same result will be yielded even when the technical idea 
of the present embodiment is applied to any one of the 
eighth, ninth, twentieth, and twenty-first embodiments. 

10 (Twenty-third Embodiment) 

An EMI analysis method according to a tenth 
embodiment of the present invention will now be described. 
In the fourteenth through twenty-second embodiments, an 
instantaneously-changing current of each cell is 

15 calculated from event information concerning each cell, 
the information being output from a logic simulator, by 
means of characterizing a peak current of a cell according 
to the type of cell and the amount of instantaneously- 
changing current for each event. The thus-calculated 

20 instantaneously-changing current is modeled as a 

rectangular waveform such that a peak value becomes equal 
to a peak current of the cell , thereby determining change 
in electric current. However, so long as the waveform of 
an instantaneous electric current is modeled as a 

25 triangular wave, a more realistic change in electric 
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current can be determined. 

By way of" example, there will now be described 
an EMI analysis method comprising the steps of: 

modeling an instantaneously-changing electric 
5 current as a triangular wave whose height becomes equal 
to a characterized peak current, the width of the 
triangular wave being calculated such that the 
instantaneously-changing current becomes equal to the 
width of the triangular wave, to thereby calculate a more 
10 realistic change in electric current; and 

subjecting the result of calculation of an electric 
waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
of interest. 

15 

FIG. 96 shows an instantaneous electric current and 
a modeled triangular wave according to the present 
embodiment. In the drawing, "i" denotes an 
instantaneously-changing current; and w t" denotes a 

20 current processing time. Further, "h" denotes the height 
of a triangular wave serving as a model of an 
instantaneously-changing current. The height of a 
triangular waveform is extracted from a library in which 
the heights of triangular waveforms of cells are 

25 characterized. Reference symbol "w" denotes the width of 
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a triangular waveform and is calculated on the basis of 
"i" and "h." An electric current at time t+x is expressed 
by Eq.5 described previously. Further, an electric 
current at time t-x is also expressed by Eq.5. 
5 The overall flow of processing of the EMI analysis 

method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 97 is a detailed block diagram showing the 

10 current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell , the 

15 information being output from the logic simulator 4 . 

Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an ins tantaneously- 

20 changing current (hereinafter referred to as an 

"instantaneously-changing current cal culation resul t 
12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 
instantaneously-changing current calculation result 12. 

25 As can be seen from FIG. 3, according to the 
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instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
5 identical with the current wave calculation operation of 
the first embodiment. Reference numeral 130 designates 
a peak current library including peak currents of cells. 
Reference numeral 160 denotes a triangular waveform height 
calculation section for calculating the height of a 
^ 10 triangular wave from the peak current library 130. 

Reference numeral 161 designates a triangular waveform 
width calculation section for calculating the width of a 
triangular waveform from the height extracted by the 
triangular waveform height calculation section 160 and the 
15 instantaneously-changing current, such that the area of 
a triangular wave becomes equal to an instantaneously- 
changing current. Reference numeral 162 denotes a 
triangular-waveform-shaping section for modeling the 
instantaneously-changing current result 12 as a 
20 triangular wave having a constant width, the height of the 
triangular waveform corresponding to the peak value of the 
triangular waveform calculated by the rectangular 
waveform height calculation section 160, and the width of 
the rectangular waveform corresponding to the width 
25 calculated by the rectangular waveform width calculation 
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section 161. 

FIG. 98 is a flowchart of processing pertaining to 
the triangular-waveform-shaping section 162, including 
processing pertaining to the triangular waveform height 
5 calculation section 160 and processing pertaining to the 
triangular waveform width calculation section 161 . First, 
the peak current library 130 is read (step 1160) . The 
following processing is iterated on a per-event basis 
until calculation of a current waveform is completed (step 

10 1161) . From the peak current library 130, a peak current 
of an instance which is an object of EMI analysis is 
extracted, and the thus-extracted peak current is taken 
as the height "h" of a triangular waveform (step 1162) . 
This processing corresponds to processing of the 

15 triangular waveform height calculation section 160. 

Provided that "w" is taken as the width of a triangular 
waveform, the area of a triangular wave is Wxh/2 and is 
equal to the instantaneously-changing current of a cell 
of the event. Therefore, the width "w" of a triangular 

20 waveform is calculated from 2xi/h (step 1163) . Reference 
symbol "i" denotes an instantaneously-changing current of 
a cell of an event, the event being an object of EMI analysis . 
This processing corresponds to processing pertaining to 
the triangular waveform width calculation section 161. 

25 Until variable "x" changes from 0 to W/2 , h(c, i) expressed 
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by Eq. Sis added to I (t+x) and I ( t-x) . Further , At is added 
to variable "x" (steps 1164 and 1165). Here, I (t+x) 
denotes total electric current flowing through all the 
cells at time t+x, and I (t-x) denotes total electric 
5 current flowing through all the cells at time t-x. This 
processing corresponds to processing pertaining to the 
triangular-waveform-shaping section 162 . 



The current wave calculation result 8 can be obtained 
10 by means of the flow of a triangular waveform height 
calculation operation, a triangular waveform width 
calculation operation, and a triangular-waveform-shaping 
operation. The FFT processing section 9 subjects the 
current waveform calculation result 8 to an FFT operation, 
15 to thereby determine the frequency characteristic of an 
EMI component. As a result, there can be obtained an EMI 
analysis result 10 . 

(Twenty- fourth Embodiment) 

20 An EMI analysis method according to a twenty-fourth 

embodiment of the present invention will now be described. 
In the fourteenth through twenty-second embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell , 

25 the information being output from a logic simulator, by 
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means of characterizing a peak current of a cell according 
to the type of cell and the amount of instantaneously- 
changing current for each event. The thus -calculated 
instantaneously-changing current is modeled as a 
5 rectangular waveform such that a peak value becomes equal 
to a peak current of the cell, thereby determining change 
in electric current. However, so long as the waveform of 
an instantaneous electric current is modeled as a Gauss 
function wave, a more realistic change in electric current 

10 can be determined. By way of example, there will now be 
described an EMI analysis method comprising the steps of: 

modeling an instantaneously-changing electric 
current as a mul ti -order-function waveform whose height 
becomes equal to a characterized peak current, the width 

15 of the multi-order-function wave being calculated such 
that the instantaneously-changing current becomes equal 
thereto, to thereby calculate a more realistic change in 
electric current; and 

subjecting the result of calculation of an electric 

20 waveform to FFT processing, thereby determining the 

frequency characteristic of an EMI component of a circuit 
of interest. 

FIG. 99 shows an instantaneous electric current and 
25 a modeled multi-order-function according to the present 
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embodiment. In the drawing, "i" denotes an 
instantaneously-changing current; and *t" denotes a 
current processing time. Further, "h" denotes the height 
of a multi-order-function wave serving as a model of an 
5 instantaneously-changing current. The height of a 

triangular waveform is extracted from a library in which 
the heights of multi-order-function waveforms of cells are 
characterized. Reference symbol "1" denotes the width of 
a multi -order-function wave and is calculated on the basis 

10 of "h," which has been calculated previously, and x ^i . " 
Provided that the width of the multi-order-function 
waveform is taken as 2xw, an electric current at time t+x 
is expressed by Eq.6. The area of the multi-order- 
function waveform is expressed by Eq.8. The area of the 

15 multi-order-function waveform is equal to the value of 
"i." Here, "a" is expressed by Eq.9. Since the height 
"h" of a multi -order-function waveform is determined by 
Eq.6 when x=0 , the height "h" is expressed by Eq.7. From 
Eq . 7 and Eq.9, the width "w" is expressed by Eq . 14 provided 

20 below. 

w = ±*ZL m (14) 
\6xh 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
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by the same block diagram as that employed in the first 
embodiment (see FIG. 1). 

FIG. 100 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
5 present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
calculating a momentary change in the current of each cell 
from the event information 5 concerning each cell, the 
information being output from the logic simulator 4. 

10 Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
conventional logic simulation. Reference numeral 12 
designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 

15 ^instantaneously-changing current calculation result 

12") determined by the instantaneously-changing current 
calculation section 11. FIG. 3 shows an example 

instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 

20 instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 5 function. 
Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 

25 the first embodiment. 
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Reference numeral 130 designates a peak current 
library including peak currents of* cells which are 
characterized according to the type of cell. Reference 
numeral 170 denotes a multi -order-function waveform 
height calculation section for calculating the height of 
a multi-order-function waveform from the peak current 
library 130. Reference numeral 171 designates a 
multi-order-function waveform width calculation section 
for calculating the width of a multi-order-function 
waveform from the height extracted by the multi-order- 
function waveform height calculation section 170 and the 
instantaneously-changing current, such that the area of 
a multi-order-function waveform becomes equal to an 
instantaneously-changing current. Reference numeral 172 
denotes a multi -order-function waveform-shaping section 
for modeling the instantaneously-changing current result 
12 into a multi -order-function waveform having a constant 
width, the height of the multi-order-function waveform 
corresponding to the peak value of the multi-order- 
function waveform calculated by the multi-order-function 
waveform height calculation section 170, and the width of 
the multi -order-function waveform corresponding to the 
width calculated by the mul ti -order-function waveform 
width calculation section 171. 
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FIG. 101 is a flowchart of processing pertaining to 
the mul ti -order-function waveform-shaping section 172, 
including processing pertaining to the mul ti-order- 
5 function waveform height calculation section 170 and 
processing pertaining to the mul ti -order-function 
waveform width calculation section 171. First, the peak 
current library 130 is read (step 1170) . The following 
J?% processing is iterated on a per-event basis until 

10 calculation of a current waveform is completed (step 1171) . 
Hi From the peak current library 130, a peak current of an 

nj instance which is an object of EMI analysis is extracted, 

O and the thus-extracted peak current is taken as the height 

D "h" of a mul ti -order-function waveform (step 1172) . This 

Q 15 processing corresponds to processing of the multi- 

order-function waveform height calculation section 170. 
The width M w" of a mul ti -order-function waveform is 
calculated by Eq.14 (step 1173). Reference symbol "i" 
denotes an instantaneously-changing current of a cell of 
20 an event, the event being an object of EMI analysis. This 
processing corresponds to processing pertaining to the 
multi-order-function waveform width calculation section 
171. Until variable "x" changes from -w to w, "y" 
expressed by Eq. 6 is added to I (t+x) and I (t-x) . Further, 
25 At is added to variable "x" (steps 1174 and 1175) . Here, 
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I(t+x) denotes total electric current flowing through all 
the cells at time t+x, and I(t-x) denotes total electric 
current flowing through all the cells at time t-x. This 
processing corresponds to processing pertaining to the 
multi-order-function waveform-shaping section 172. 
[0255] 

The current wave calculation result 8 can be obtained 
by means of the flow of a multi-order-function waveform 
height calculation operation, a multi -order-function 
waveform width calculation operation, and a multi- 
order-function waveform-shaping operation. The FFT 
processing section 9 subjects the current waveform 
calculation result 8 to an FFT operation, to thereby 
determine the frequency characteri stic of an EMI component. 
As a result, there can be obtained an EMI analysis result 
10. 

(Twenty-fifth Embodiment) 

An EMI analysis method according to a twenty-fifth 
embodiment of the present invention will now be described. 
In the fourteenth through twenty-second embodiments, an 
instantaneously-changing current of each cell is 
calculated from event information concerning each cell, 
the information being output from a logic simulator, by 
means of characterizing a peak current of a cell according 
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to the type of cell and the amount of instantaneously- 
changing current for each event. The thus-calculated 
instantaneously-changing current is modeled as a 
rectangular waveform such that a peak value becomes equal 
5 to a peak current of the cell, thereby determining change 
in electric current. However, so long as the waveform of 
an instantaneous electric current is modeled as a Gauss 
function waveform, a more realistic change in electric 
current can be determined. 

10 By way of example, there will now be described an 

EMI analysis method comprising the steps of: 

modeling, for each type of cell, an 
instantaneously-changing electric current as a Gauss 
function waveform whose height and width are characterized, 

15 to thereby calculate a more realistic change in electric 
current; and 

subjecting the result of calculation of an electric 
waveform to FFT processing, thereby determining the 
frequency characteristic of an EMI component of a circuit 
20 of interest. 

FIG. 50 shows an instantaneous electric current and 
a modeled Gauss function according to the present 
embodiment. In the drawing, "i" denotes an 
25 instantaneously-changing current; "t" denotes a current 
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processing time; and *h" denotes the height of a Gauss 
function wave serving as a model. The height "h" is equal 
to a peak current characterized according to the type of 
cell. Reference symbol x> w" denotes a width of a Gauss 
5 function serving a model, and the width M w" corresponds 
to a value characterized according to the type of cell. 
In the present embodiment, a Gauss function serving as a 
model is expressed as Eq.ll. This function is called a 

JQ Gaussian function. Let x=0 , Eq.ll provides "h," which 

on 

M 10 represents the height of a Gaussian function waveform. 

Ul From Eq . 11 , let x=3xW , there is derived a value of 0 . OOOlxh 

^ (i.e., 0.01% of the height of a Gaussian function wave) . 

H This value is taken as a limit on the width of a Gaussian 

~ function serving as a model, and the width of the Gaussian 

™ 15 function is taken as xx 6xw." 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 
by the same block diagram as that employed in the first 
20 embodiment (see FIG. 1) . 

FIG. 102 is a detailed block diagram showing the 
current waveform calculation section 7 according to the 
present embodiment. Reference numeral 11 designates an 
instantaneously-changing current calculation section for 
25 calculating a momentary change in the current of each cell 
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from the event information 5 concerning each cell, the 
information being output from the logic simulator 4 . 
Calculation of momentary change in current is technically 
analogous to calculation of power consumption through a 
5 conventional logic simulation. Reference numeral 12 

designates a calculation result of an instantaneously- 
changing current (hereinafter referred to as an 
^instantaneously-changing current calculation result 
12") determined by the instantaneously-changing current 

10 calculation section 11. FIG. 3 shows an example 

instantaneously-changing current calculation result 12. 
As can be seen from FIG. 3, according to the 
instantaneously-changing current calculation result, the 
waveform of an electric current assumes a 8 function. 

15 Therefore, an EMI component cannot be analyzed through an 
FFT operation. Thus far, the processing has been 
identical with the current wave calculation operation of 
the first embodiment. A Gaussian function waveform 
height calculation section 180 extracts a peak value of 

20 a Gaussian function serving as a model, from the peak 
current library 130 including peak currents of cells 
characterized according to the type of cell. Reference 
numeral 181 designates a wave width library includes wave 
widths of cells characterized according to the type of cell . 

25 Reference numeral 182 designates a Gaussian function 
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waveform width calculation section for extracting, from 
the wave width library 181, the width of a Gaussian 
function waveform serving as a model . 

5 Reference numeral 183 designates a Gaussian function 

waveform- shaping section for modeling the 
instantaneously-changing current result 12 as a Gaussian 
function wave having a peak value calculated by the 
Gaussian function waveform height calculation section 180 
10 and a width calculated by the Gaussian function waveform 
width calculation section 182. 

FIG. 103 shows the flowchart of processing 
pertaining to the Gaussian function waveform-shaping 

15 section 183, including processing pertaining to the 

Gaussian function waveform height calculation section 180 
and processing pertaining to the Gaussian function 
waveform width calculation section 182. First, the peak 
current library 130 is read (step 1180) , and the wave width 

20 library 181 is read (step 1181) . 

The processing provided below is iterated on a 
per-event basis until calculation of a current waveform 
is completed (step 1182) . From the peak current library 
25 130, a peak current of an instance which is an object of 
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EMI analysis is extracted, and the thus-extracted peak 
current is taken as the height u h" of a Gaussian function 
waveform (step 1183) . This processing corresponds to 
processing of the Gaussian function waveform height 
5 calculation section 180. The width "w" of a Gaussian 
function waveform which is an object of EMI analysis is 
extracted from the wave width library 181, and the 
thus-extracted width is taken as the width M w" of a 

Q 

y3 Gaussian function waveform (step 1184) . This processing 

10 corresponds to processing of the Gaussian function 

ITr i 

in waveform width calculation section 182. Until variable 

l ^ "x" changes from -3xw to 3xw, "y" expressed by Eq.ll is 

™ added to I(t+x). Further, At is added to variable M x" 

(steps 1185 and 1186). Here, I(t+x) denotes total 
~ 15 electric current flowing through all the cells at time t+x . 

This processing corresponds to processing pertaining to 
the Gaussian function waveform-shaping section 183. 

The current wave calculation result 8 can be obtained 
20 by means of the flow of a Gaussian function waveform height 
calculation operation, a Gaussian function waveform width 
calculation operation, and a Gaussian function 
waveform-shaping operation. The FFT processing section 
9 subjects the current waveform calculation result 8 to 
25 an FFT operation, to thereby determine the frequency 
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characteristic of an EMI component. As a result, there 
can be obtained an EMI analysis result 10. 

(Twenty- sixth Embodiment) 
5 An EMI analysis method according to a twenty-sixth 

embodiment of the present invention will now be described. 

The EMI analysis method according to the present 
embodiment is characterized by modeling an 
instantaneously-changing current while the information 

10 concerning a current waveform of each circuit element (or 
cell) which is described in connection with the fourteenth 
through twenty-second embodiments is separated into an 
element dependent on a short circuit current (also called 
"short circuit-current-dependent element") and another 

15 element dependent on a charge current (also called 

"charge-curren t-dependent element") . In contrast with 
the thirteenth embodiment, the present embodiment effects 
modeling of each of the circuit elements (cells) , thereby 
enabling formation of a more accurate model of electric 

20 current. 

For example, in a case where the circuit connection 
information shown in FIG. 110 and the pattern of signal 
change shown in FIG. 12 are imparted to a circuit element, 
25 an electric current shown in FIG. 104 flows through the 
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circuit element. In such a case, the conventional EMI 
analysis method enables formation of only a waveform such 
as that shown in FIG. 14 . Even in the case of the 
thirteenth embodiment, there is produced only a waveform 
such as that shown in FIG. 58. In contrast, the EMI 
analysis method according to the present embodiment 
enables formation of a more correct electric current . FIG . 
72 shows the configuration of an EMI analysis method 
according to the present embodiment, the method using the 
value of a peak current and the amount of electric current 
described in connection with the fourteenth embodiment. 

A system to be used with the EMI analysis method shown 
in FIG. 72 comprises circuit connection information 
storage means 9001; pattern-of -signal -change storage 
9002; element current storage means 9003 for storing the 
amount of electric current flowing in a circuit element; 
total current storage means 9004 for storing a total amount 
of electric current flowing through the circuit element; 
and total current calculation means 9005 for calculating 
the total amount of electric current flowing through a 
circuit element. 

Next will be described individual elements 
constituting the EMI analysis method shown in FIG. 72 and 
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procedures for calculating the total amount of electric 
current shown in FIG. 106, from the circuit connection 
information shown in FIG . 110 , the pattern of signal change 
shown in FIG. 12, and the amount of electric current 
flowing through a circuit element shown in FIG. 105. 

The circuit connection information storage means 
9001 stores information concerning the connection of a 
circuit to be analyzed in the same manner as described in 
connection with the background art. Circuit connection 
information, such as that shown in FIG. 110, is stored in 
advance in the circuit connection information storage 
means 9001. 

The pattern-of -signal -change storage 9002 stores 
the pattern of signal change to be applied to the circuit 
connection information, as in the case of a conventional 
pattern of signal change described in connection with the 
background art. As shown in FIG. Ill, a pattern of change 
in a signal input to an input port CK and a pattern of change 
in a signal input to an input port D, such as those shown 
in FIG. 110, are stored in the pattern-of -signal -change 
storage 9002 beforehand. 

The element current storage means 9003 stores^.. for 
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each circuit element and in a form such as that shown in 
FIG. 105, the total amount of electric current flowing in 
a circuit element (cell) , a peak value of the electric 
current, the total amount of charge current, and a peak 
5 current. 

The total current calculation means 9005 operates 
according to a flowchart, such as that shown in FIG. 107. 
According to the flowchart relating to a method of 

10 calculating an electric current shown in FIG. 107, 

processing pertaining to step 1504 of the conventional 
current calculation method is replaced by processing of 
calculating a rectangular waveform of the short circuit 
current from a peak value and the total amount of electric 

15 current, calculating a rectangular waveform of the charge 
current from a peak value and the total amount of electric 
current, and merging the thus -calculated rectangular 
waveforms into a single rectangular waveform (step 12504) . 
In step 12501, there is read the circuit connection 

20 information which is stored in the circuit connection 
information storage means 9001 and shown in FIG. 110. 



In step 12502, ther 
change which is shown in 
25 pattern -of -signal -change 



e is read a pattern of signal 
FIG. Ill and is stored in the 
storage 9002. In step 12503, 

181 



there is read the element current information which is 
shown in FIG. 105 and stored in the element current storage 
means 9003. In step 12504, the pattern of signal shown 
in FIG. Ill is imparted to the circuit connection 
5 information shown in FIG. 110, thereby enabling 

propagation of a signal. In a case where a change arises 
in a signal appearing at an external input terminal of a 
circuit element, which element is described in the 
information concerning the amount of electric current 

10 flowing in an element, a rectangular waveform of an 

electric current is produced by means of averaging a short 
circuit current component of the electric current by a peak 
value of the short circuit current component, and a 
rectangular waveform of an electric current is produced 

15 by means of averaging a charge current component of the 
electric current by a peak value of the charge current 
component. The thus-formed rectangular waveforms are 
added to the information concerning the total mount of 
electric current obtained at a time when the change has 

20 arisen in the signal appearing at the external input 

terminal. After completion of propagation of a signal, 
in step 12505 the thus-calculated information concerning 
the total amount of electric current obtained at each time 
shown in FIG. 106 is stored in the total current storage 

25 means 9004. As a result, there is produced a current 
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waveform, such as that shown in FIG. 106, which is more 
realistic than that which is obtained by the conventional 
current calculation method and shown in FIG. 14. 



5 (Twenty-seventh Embodiment) 

Next will be described an EMI analysis method 
according to a twenty-seventh embodiment of the present 
invention. The tenth embodiment has described that a more 
realistic change in current can be materialized by means 

10 of calculating an instantaneously-changing current of 
each cell from event information concerning the cell, the 
information being output from a logic simulator, and 
modeling the thus -calculated instantaneously-changing 
current as a triangular wave. The EMI analysis method 

15 according to the present embodiment is characterized by 
calculating a current waveform in consideration of output 
slew information; that is, the time which elapses from the 
time a change arises in a signal appearing at an output 
terminal of a cell until the change has been canceled. 

20 By way of example, there will now be described an 

EMI analysis method comprising the steps of: 

modeling an instantaneously-changing current as a 
triangular waveform by means of calculating the width of 
a triangular waveform for each event information in 

25 consideration of slew information concerning an output 
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terminal of a cell and the height of the triangular 
waveform on the basis of the thus-calculated width such 
that the area of the triangular waveform becomes equal to 
the amount of current of each event; 
5 calculating a more realistic change in electric 

current by means of using an existing library in its 
present form; and 

subjecting the result of calculation of an electric 
y3 waveform to FFT processing, thereby determining the 

M 8 10 frequency characteristic of an EMI component of a circuit 

in! 5 

of interest. 



FIG. 117 shows a logic value (see FIG. 117A) and a 
modeled triangular waveform (see FIG. 117B) , which are to 

15 be used for executing the EMI analysis method according 
to the present embodiment. Reference symbol "t" denotes 
a current processing time, and ^W" denotes the width of 
a triangular waveform serving as a model of an 
instantaneously-changing current. The width "W" is 

20 determined by means of calculating an instantaneously- 
changing current of each cell from event information of 
each cell, in consideration of output slew information. 
Further, M h" denotes the height of a triangular waveform 
and is calculated from "i" and «W." Here, h=2xI/W. 

25 FIG. 118 shows the configuration of a system for 
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effecting the EMI analysis method according to the present 
invention . The system used for effecting the EMI analysis 
method shown in FIG. 118 comprises pattern-of -signal - 
change storage means 1171; element current storage means 
5 1172 for storing the amount of electric current flowing 
in a circuit element; output slew information storage 
means 1173; and current waveform calculation means 1175 
for calculating a current waveform on the basis of data 
output from the pattern-of -signal-change storage means 

10 1171, the element current storage means 1172, the output 
slew information storage means 1173, and circuit 
connection information storage means 1174 . The system 
outputs a calculation result of a current waveform 1176 
(hereinafter referred to simply as a "current waveform 

15 calculation result 1176") . 

FIG. 119 is a flowchart of a triangular- 
waveform-shaping operation, including a triangular 
waveform height calculation operation and a triangular 

20 waveform width calculation operation. The following 

processing is iterated on a per-event basis after reading 
of an output slew (step 1190) until calculation of a 
current waveform is completed (step 1191) . Abase of the 
triangular waveform is extracted as an output slew while 

25 a point at which a change arises in an event is taken as 
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a center (step 1192). An isosceles triangle whose area 
is equal to the amount of" electric current flowing through 
the event is taken as a current waveform of the event. The 
area of a triangular wave is Wxh/2 and assumes a value M i." 
5 The height "h" of a triangular wave is calculated from 
2xi/W (step 1193) . Reference symbol "i" denotes an 
instantaneously-changing current of a cell of an event 
which is an object of EMI analysis. This processing 
corresponds to processing pertaining to the triangular 

10 waveform height calculation section 100 according to the 
tenth embodiment. Until variable "x" changes from 0 to 
W/2, (2xh/W)*x+h is added to I (t+x) and I (t-x) . Further, 
At is added to variable "x" (steps 1194, 1195, and 1196) . 
Here, I (t+x) denotes total electric current flowing 

15 through all the cells at time t+x, and I (t-x) denotes total 
electric current flowing through all the cells at time t-x. 
This processing corresponds to processing pertaining to 
the triangular-waveform-shaping section 101. 

20 The current wave calculation result 76 can be 

obtained by means of the flow of a triangular waveform 
height calculation operation and a triangular- 
waveform-shaping operation. The FFT processing section 
9 subjects the current wave calculation result 76 to an 

25 FFT operation, to thereby determine the frequency 
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characteristic of an EMI component. As a result, there 
can be obtained an EMI analysis result 10. 

The EMI analysis method enables estimation of an 
5 electric current with sufficient accuracy from only the 
information which is currently used by a power consumption 
analysis tool and used for verifying a delay. 

(Twenty- eighth Embodiment) 

10 Next will be described an EMI analysis method 

according to a twenty-eighth embodiment of the present 
invention. The EMI analysis method according to the 
present embodiment is characterized by means of 
calculating a current waveform, in consideration of 

15 whether an event is a rise or a fall and in addition to 
the feature of the EMI analysis method according to the 
twenty-seventh embodiment. FIG. 120 shows a logic value 
(see FIG. 120A) and a modeled triangular waveform (see FIG. 
120B) , which are to be used for effecting an EMI analysis 

20 method according to the present embodiment. Reference 
symbol M t" denotes a current processing time, and XX W" 
denotes the width of a triangular waveform serving as a 
model of an instantaneously-changing current. The width 
is determined by means of calculating an 

25 instantaneously-changing current of each cell from event 
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information of each cell, in consideration of output slew 
information. There is employed rise-dependent 
coefficient /fall -dependent coefficient storage means 
1207, and current waveform calculation means 1205 
5 calculates a current waveform in accordance with whether 
an event is a rise or a fall. 

The overall flow of processing of the EMI analysis 
method according to the present embodiment is represented 

10 by the same block diagram as that employed in the first 
embodiment (see FIG. 1) . 

FIG. 121 shows the configuration of a system for 
effecting the EMI analysis method according to the present 
invention . The system used for effecting the EMI analysis 

15 method shown in FIG. 121 comprises pa ttern-of -signal - 
change storage means 1201; element current storage means 
1202 for storing the amount of electric current flowing 
in a circuit element; output slew information storage 
means 1203; and current waveform calculation means 1205 

20 for calculating a current waveform on the basis of data 
output from the pattern-of -signal-change storage means 
1201, the element current storage means 1202, the output 
slew information storage means 1203, and circuit 
connection information storage means 1204 . The system 

25 outputs a calculation result of a current waveform 1206 
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(hereinafter referred to simply as a "current waveform 
calculation result 1206") . 

In the present embodiment, a current waveform is 
5 calculated in consideration of whether an event is a rise 
or a fall, by reference to the rise/fall -dependent 
coefficient storage means 1207, in addition to the 
operations of the EMI analysis method described in 
connection with the twenty-seventh embodiment. A more 
10 realistic change in current is calculated by use of an 
existing library in its present form, and by subjecting 
a calculation result of a current wave to FFT processing, 
to thereby determine the frequency characteristic of an 
EMI component and enable EMI analysis. 



according to a twenty-ninth embodiment of the present 
invention. According to the present method, there is 

20 estimated the amount of electric current on the basis of 
event information, the number of cells determined 
according to a SPICE netlist provided in a cell, an output 
load capacitance, input slew information, the 
relationship between an input slew and the amount of 

25 current, and a time at which a change has arisen in an input 
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( Twenty-ninth Embodiment) 



Next will be described an EMI analysis method 
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signal, or a delay time. Information concerning the 
amount of electric current which flows in a single event 
is calculated from an output load capacitance and a table 
representing the relationship between an input slew and 
5 the amount of electric current. FIG. 122 shows a logic 
value (see FIG. 122A) and a modeled triangular waveform 
(see FIG. 122B) , which are to be used for effecting an EMI 
analysis method according to the present embodiment. 
Reference symbol *t d " designates a delay time; w t ff denotes 

10 a current processing time; and "W" denotes the width of 
a triangular waveform serving as a model of an 
instantaneously-changing current. The width "W" is 
determined by means of calculating an instantaneously- 
changing current of each cell from event information of 

15 each cell, on the basis of output slew information taking 
into consideration a delay time td caused by multistage 
connection. Here, there is calculated the amount of 
current independent of an output load (i.e., the amount 
of non-dependent current) beforehand, and the thus- 

20 calculated amount of current is stored 0 . A base of the 
triangular waveform is taken as an output slew while a 
point at which a change arises in an output is taken as 
a center. Further, there are employed, as output loads, 
an angular waveform whose area is equal to the amount of 

25 non-dependent current or to the amount of current flowing 
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electric current-capacitance is 0 and a triangular 
waveforms equal in number to the number of cells within 
the delay time. As shown in FIG. 123, the EMI analysis 
5 method according to the present embodiment is provided 
with a load capacitance/slew-dependent element current 
information table 122 6; element current calculation means 
1227 for calculating the amount of electric current 
flowing through a cell from the load capaci tance/slew- 

10 dependent element current information table 122 6; and 
non-load element current storage means 1228 for storing 
the amount of electric current flowing through a non-load 
element calculated by the element current calculation 
means 1227. By use of these means, a current waveform is 

15 calculated. The systemwhich is shown in FIG. 123 and used 
for effecting the EMI analysis method according to the 
present invention, comprises the load 

capaci tance/slew-dependent element current information 
table 1226; the element current calculation means 1227 for 

20 calculating the amount of electric current flowing through 
a cell from the load capacitance/slew-dependent element 
current information table 1226; the non-load element 
current storage means 1228 for storing the amount of 
electric current flowing through a non-load element 

25 calculated by the element current calculation means 1227; 
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pattern-of -signal -change storage means 1221; element 
current storage means 1222 for storing the amount of 
electric current flowing in a circuit element at the time 
of an event, the electric current being calculated by the 
5 element current calculation means 1227; output slew 
information storage means 1223; and current waveform 
calculation means 1225 for calculating a current waveform 
on the basis of data output from the pattern-of- 
signal-change storage means 1221, the element current 

10 storage means 1222, the output slew information storage 
means 1223, circuit connection information storage means 
1224, number-of -stages storage means 1229 for storing the 
number of stages in a composite cell, and the non-load 
element current storage means 1228 . The system outputs 

15 a calculation result of a current waveform 1206 

(hereinafter referred to simply as a "current waveform 
calculation result 1206") . 



In the twenty-ninth embodiment, there are stored a 
20 plurality of triangular waveforms which are equal in 

number to cells, such that the total area of a triangular 
waveform becomes equal to the amount of electric current 
obtained when a capacitance assumes a value of 0, during 
a period of time from the time when a change has arisen 
25 in an input signal until a change arises in an event. A 
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more realistic change in current is calculated from an 
existing library in its present form, and the calculation 
result of current waveform is subjected to FFT processing, 
thereby determining the frequency characteristic of an EMI 
5 component arising in a circuit to be analyzed and enabling 
EMI analysis. 



The present invention evaluates electromagnetic 
10 interference arising in an LSI through a simulation, by 
means of accurately analyzing a change in a power current, 
which change can be said to primarily account for 
electromagnetic interference. In order to solve the 
problem, a highly-accurate power current analysis 
15 function is added to a gate-level simulation, thereby 
preventing consumption of much time, which would otherwise 
be required by the conventional EMI analysis method using 
transistor-level current analysis. 
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